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survey,  conducting  the  survey  and  reporting  on  the  survey  were  included  In  the 
study. 

The  site  modeling  program  models  the  manufacturing  processes  which  contri 
bute  pollutants  to  the  system,  models  the  sewer  system,  and  models  the  treat- 
ment system  including  acid  or  caustic  neutralization,  settling  ponds,  and 
domestic  treatment.  The  inputs  to  the  model  are  the  production  levels  of  the 
manufacturing  .processes  and  the  outputs  are  the  predicted  pollutant  measuremen 
values  at  each  possible  measure  point  in  the  system. 

The  resource  matching  program  accepts  data  defining  proposed  measurements 
and  matches  these  against  the  available  time,  manpower,  and  equipment.  The 
output  lists  the  pollutant  to  be  measured  at  each  measure  point,  the  total 
commitment  of  time  for  each  analyst  and  for  each  piece  of  equipment.  Note  is 
made  of  any  overcommitment  of  manpower  or  equipment. 

The  model  refinement  or  updating  program  accepts  measurements  taken  during 
a preliminary  survey  or  during  a regular  survey  and  computes  suggested  new 
parameters  for  the  process  models. 

The  indicator  model  program  evaluates  the  performance  of  sanitarv  treat- 
ment facilities. 

The  program  uses  design  data,  data  from  the  operating  log  and/or  data 
generated  during  the  survey  and  computes  key  operational  characteristics.  Com 
paring  these  with  desirable  values  as  cited  in  design  books  and  manuals  will 
give  the  survey  planner  insight  into  the  operation  of  the  system  and  suggest 
the  need  for  more  survey  measurements  or  the  need  for  changes  in  operation. 

A system  was  developed  for  automatic  instrumentation  of  pH,  conductivity, 
and  other  parameters  which  use  strip  charp  recordings.  Interface  hardware  was 
selected  and  purchased  and  Interface  software  was  developed  for  direct  connec- 
tion to  a digital  computer. 

A data  handling  system  was  developed  for  use  during  and  after  the  survev. 

A PDP8-OS/8  and  peripheral  equipment  was  purchased.  Software  was  developed  to 
perform  data  handling  functions  and  to  direct  the  user  in  application  of  the 
software.  The  program  accepts  raw  data  from  the  analytical  chemist  and  per- 
forms data  conversions,  transcriptions,  and  data  logging  functions.  Output 
is  available  in  several  forms  as  may  be  needed  for  various  reports  during  and 
at  the  end  of  the  survey. 

Recoimnendations  are:  the  survey  planner  should  obtain  sufficient  dat.i  in 

a prollminary  survey  to  model  and  analyse  the  site;  measurements  should  he  .aut  •- 
mated  to  the  maximum  extent  possible;  data  handling  should  be  delegated  to  the 
computer  when  the  operations  are  well  defined  and  repetitive.  The  programs, 
software  and  hardware  included  here  will  assist  the  survey  planner  in  followin 
these  recommendations  and  design  a more  effective  survey. 
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complete  so  that  an  experienced  programmer  should  be  able  to  add  to, 
delete  or  alter  parts  of  this  program. 

For  each  subroutine  the  variable  names  and  variable  notations  are 
defined.  The  function  of  each  statement  or  small  group  of  statements 
is  explained.  Flow  diagrams  are  included  to  depict  the  logic  relation- 
ship among  statements  and/or  subroutines.  . 


CHANGING  SYSTEM  PLANNING  PROGRAM  DIMENSIONS 


The  program  and  subroutines  have  been  written  so  that  only  cards 
in  the  MAIN  program  need  to  be  changed  to  redimension  the  system.  The 
quantities  that  can  be  varied  are;  the  number  of  pollutant  parameters, 
the  number  of  sources  (processes),  the  number  of  branches,  the  number 
of  pieces  of  equipment. 

The  specific  items  in  the  MAIN  program  which  must  be  changed  to 
redimension  the  system  are  given  in  the  following  paragraphs: 

To  handle  more  than  25  sources  (where  i = number  of  sources) 
in  the  MAIN  program, 

change  MSORS  = 25  to  MSORS  = i and  in 

the  DIMENSION  statement  the  arrays  to  modify  are: 


XNAME  (25) 

YNAME  (25) 

FLOW  (25) 

SPLIST  (25,25)  CHANGE  TO 

CAP  (25) 

A (25,25) 

C (25,25) 

P (25,25) 

SUM  (25) 

SUMA  (25) 


XNAME  (i) 
YNAME  (i) 

FLOW  (i) 
SPLIST  (i,25) 
CAP  (i) 

A (i,25) 

C (i,25) 

P (i,25) 

SUM  (i) 

SUMA  (i) 


To  handle  more  than  25  branches  (where  n = number  of  branches) 


in  the  program. 


change  MBRNC  = 25  to  MBRNC  = N 


and  in  the  DIMENSION  Statement  the  arrays  to  modify  are; 


A (25,25) 

Y (25,25) 

X (25,25) 

EFF  (25) 

NBRNCH  (25,2)  CHANGE  TO 

I BN  (51,26) 

PC  (25,25,25) 

AMAR  (25,27) 


A (25, n) 

Y (n,  25) 

X (n,25) 

EFF  (n) 

NBRNCH  (n,2) 
IBN  (5i,n+l) 
PC  (25,25,n) 
AMAR  (25,n+2) 
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NALOW  (25,25) 
BRN  (25) 

BRANCH  (25,25) 
NROUT  (25) 
NFLOW  (25) 


NALOW  (25, n) 
BRN  (n) 

BRANCH  (n,25) 
NROUT  (n) 
NFLOW  (n) 


To  handle  more  than  25  parameters  (m  = number  of  parameters), 
in  the  MAIN  program, 

change  MPARM  = 25  to  MPARM  = m 

and  in  the  DIMENSION  statement  the  arrays  to  modify  are; 


SPLIST  (25,25) 

SPLIST  (25, m) 

c (25,25) 

C (25, m) 

Y (25,25) 

Y (25, m) 

X (25,25) 

X (25, m) 

P (25,25) 

P (25, m) 

NPLIST  (25) 

NPLIST  (m) 

POLN  (25,5) 

POLN  (m,5) 

NTEMP  (25) 

NTEMP  (m) 

SAMFRE  (25) 

SAMFRE  (m) 

NMA  (25)  rwANCF 

MENAME  (25,3,5) 

TO 

NMA  (m) 

MENAME  (m,3,5) 

PCRM  (25,3) 

PCRM  (m,3) 

PMDATA  (25,3,4,5) 

PMDATA  (m,3,4,5) 

IDO  (25,25) 

IDO  (m,25) 

PM  (25,25,3) 

PM  (m,25,3) 

AMAR  (25,27) 

AMAR  (m,27) 

PC  (25.25,25) 

PC  (m,25,25) 

NALOW  (25,25) 

NALOW  (m,25) 

EQUSED  (25,75) 

EQUSED  (m,75) 

NSET  (25,3) 

NSET  (m,3) 

FLGPT  (25) 

FLGPT  (m) 

RANK  (25) 

RANK  (m) 

BRANCH  (25,25) 

BRANCH  (m) 

ELE  (25) 

ELE  (m) 

To  handle  more  that  75 
in  the  MAIN  program. 

pieces 

of  equipment  (j  = number  of 

change  MEQ,  = 75 

to 

II 

O' 

UJ 

and  in  the  DIMENSION  statement: 

change  VNSP (75) 

to 

VNSP  (j) 

SMEQTI  (75) 

SMEQTI  (j) 

EQUSED  (25,75) 

CHANGE 

TO  EQUSED  (25. j) 
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EU  (75) 

EtiTIME  (75) 
EQNAME  (75,5) 


EU  (j) 

EQTIME  (j) 
EQNAME  (j,5) 
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It  should  be  noted  that  these  are  the  only  changes  required  on  the 
array  dimensions.  There  is  no  need  to  make  modifications  in  the  sub- 
routines. Also,  it  is  very  important  that  the  array  parameters  MSORS, 
MBRNC,  MPARM,  MEQ,  MBRPl , and  MP2  and  the  array  subscripts  correspond 
exactly.  For  example,  if  MPARM  = 31,  then  every  array  that  is  affected 
by  an  increase  in  parameter  capability  must  have  ''31"  appearing  in  the 
proper  subscript  position  as  explained  in  that  section.  Care  chould 
also  be  exercised  that  every  array  has  its  proper  subscript  changed. 

Any  omission  may  lead  to  errors  that  may  not  be  immediately  apparent. 
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MAIN  - This  program  calls  the  subroutines  LEVEL,  START,  CORRCT,  PROCES, 

TOP,  and  RM.  Also  It  sets  up  arrays  and  zeros  out  memory  space. 
Variable  Definition 

Variables  are  defined  in  each  subroutine  where  they  appear. 

Listed  here  are  only  the  variables  specifically  referred  to  in  this 
program. 

NFLAG 
MPARM 

TEST 
MSORS 
MBRNC 
MEQ 

MBRPl 
MP2 

MASTER 

EPSLON 
Statements 

1 defines  COMMON 

2-7  dimensions  all  arrays  in  the  program.  These  dimensions  must  be 

in  agreement  with  numbers  assigned  MPARM,  MSORS,  MBRNC,  MEQ, MBRPl, 
and  MP2. 

0-13  defines  size  of  MPARM,  MSORS,  MBRNC,  MEO,  MBRPl,  MP2 

14-40  zeros  out  arrays  NSET,  PCRM,  PMDATA,  SAMFRE,  NMA,  NALOW,  Y.  PC, 

PM,  AMAR,  NPLA,  SMEQTI,  EQTIME,  CONSTR,  SUMM. 


- a flag  to  determine  if  the  RM  routine  is  to  be  called 

- a number  to  dimension  arrays--the  number  of  pollutant 
parameters 

- a flag  to  indicate  insufficient  data  for  TOP 

- a number  to  dimension  arrays--the  number  of  sources 

- a number  to  dimension  arrays--the  number  of  branches 

- a number  to  dimension  arrays--the  number  of  pieces  of 

equi pment 

- the  number  of  branches  plus  1 

- the  number  of  branches  plus  2 

- the  largest  parameter  number  to  be  used  in  the  current 
problem 

- a flag  to  determine  if  CORRCT  is  to  be  called 


I 


I 


41  MASTER  is  read  in 

42  NFLAG  is  read  in 

43  EPSLON  is  read  in 

44  format  for  reading  MASTER,  NFLAG,  EPSLOf! 

45  CALL  PROCES  subroutine 

46  CALL  START  subroutine 

47  CALL  TOP  subrouting 

48  if  TEST  = 1 terminate  program  if  ALK  and  ACV  are  not  included 

in  NPLIST 

49  if  EPSLON  = 0 d-j  not  call  CORRCT 

50  Call  CORRCT 

51  if  NFLAG  = 0 call  LEVEL  and  RM 

52  Call  LEVEL 

53  Call  RM 

54  STOP 

55  END 
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SUBROUTINE  PROCES  - This  program  accesses  a library  of  subroutines 


each  of  which  defines  a process  in  terms  of  flow  at 
100^  capacity  operation  and  concentration  or  other 
measure  number  of  pollutant  parameter  from  the  process. 
It  forms  two  arrays,  one  for  flows  from  each  process 
and  another  for  concentrations  of  each  pollutant 
from  each  process.  It  also  allows  the  modification 
of  these  arrays  by  changing  individual  parameters  with- 
in the  arrays  or  by  rep  acing  entire  rows  which 
represents  an  entire  process. 

Variable  Def ini tion 


XNAME  (MSORS) 


L 


YNAME  (MSORS) 

ELE  (MPARM) 
VAL  (MPARM) 

ZNAME  (25) 


FEW  (25) 

FLOW  (MSORS) 


SPLIST  (MSORS, 
MPARM) 


CAP  (MSORS) 


- orocess  names--read  in  when  defining  list  of 
processes  in  a plant 

- process  names--read  in  wnen  performing  mod i f i ca t ions -- 
defines  row  in  which  modification  is  to  be  made 

- parameter  (element)  in  row  to  be  modified 

- value  to  which  above  parameter  (element)  is 
to  be  modified 

- process  model  names--name  of  process  models  which 

do  not  reside  in  library  which  user  wishes  to  define 

- value  which  represents  100^  flow  of  a new  process  model 

- array  which  holds  values  of  process  flows  at  full 
operating  capacity 

- array  which  holds  values  of  pollutant  parameters-- 
rows  relate  to  process  model  names--col umns  relate 
to  each  different  pollutant  name 

- operating  capacity  of  a process 
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L i nes 


I 

i 


I 1 

* 2 

1 ^ 

it-  54 

I 55-  60 

61-  62 

• 63-  68 

69 

I 70-122 

( 

j 123-124 

I 

125 

126 

127-128 

129 

130-134 

135-139 


Name  of  subroutine  PROCES 

Dimensions  the  arrays,  XNAME,  ELE,  YNAME,  VAL,  ZNAME,  FLW, 

FLOW,  SPLIST,  CAP 

Sets  the  variable  ELE  to  an  integer  variable 
Initialize  the  names  of  processes  in  the  library 
Initialize  the  FLOW  and  SPLIST  arrays  to  zero 

Reads  N - the  number  of  processes  to  be  used  in  a particular  run 
Reads  and  writes  the  names  of  each  process  to  be  used  in  a 
run  and  the  correspond i ng  operating  capacity  of  each  process. 
Resets  JM  = 0 

Selects  and  calls  the  appropriate  subroutine  for  defining 
each  process  used  in  a particular  run.  These  subroutines 
are  described  below 

Reads  NOD  - this  variable  if  "1"  indicates  that  these  are 
modifications  to  be  performed  on  the  SPLIST  array  - if  "0" 
indicates  that  there  are  no  modifications 
Wr  i tes  NOD 

Checks  value  of  NOD  - if  NOD  = "I"  program  execution  is 
transferred  to  statement  127  " if  NOD  = "0”  program  execution 
is  transferred  to  Statement  140 

Reads  NMOD  - the  number  of  modifications  to  the  SPLIST  array 
Writes  NMOD 

Reads  and  writes  YNAME,  ELE,  and  VAL  variables 
Compares  YNAME 's  to  XNAME 's.  When  YNAME  matches  XNAME  the 
value  of  the  SPLIST  element  defined  by  JM  and  ELE(NM)  is 
replaced  by  VAL(NM) 
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140-142 


143 


144-145 

146 


147-159 


160 

161 


Reads  and  writes  NPRO  - if  NPRO  = "1"  indicates  that  there 
are  processes  named  in  the  XNAME  array  which  are  not  in  the 
library  and  which  are  to  be  defined  now.  If  NPRO  = "0” 
indicated  that  there  are  no  undefined  processes  in  the  XNAME 
array. 

Checks  the  value  of  NPRO  - if  NPRO  = "I"  execution  is 
transferred  to  Statement  144.  If  NPRO  = "0"  execution 
is  transferred  to  Statement  160. 

Reads  NUMP  - the  number  of  new  processes  to  be  defined  now 
Reads  NELE  - the  number  of  pollutant  elements  to  be  placed 
in  the  SPLIST  row  pertaining  to  each  new  process  being 
defined  now. 

Does  the  following: 

1.  Reads  ZNAME  - name  of  new  process  previously  put 

in  XNAME  array  and  also  reads  FLW  which  is  flow  from 
new  process  at  1001  capacity. 

2.  Compares  ZNAME  to  XNAME  array  entries.  When  proper 
entry  found  this  process  is  defined; 

a.  FLOW  - flow  for  new  process  at  operating  capacity 

b.  SPLIST  - elements  are  read  in  corresponding  to  new 
process . 

This  process  is  repeated  for  each  new  process  until  all  new 
processes  are  defined. 

Return  statement  - returns  program  execution  to  main  program 
END  statement 
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PROCESS  LIBRARY  - Subroutines  for  defining  processes  in  library. 

One  subroutine  for  each  process  in  library. 


Line 


1 

2 

3 

Next  N 
statements 


Subroutine  name 

DIMENSIONS  FLOW,  SPLIST,  CAP 

Defines  flow  of  process  at  operating  capacity 

Flow  (JM)  = .01  * CAP  (JM)  in  % * MGD  at  100^ 

The  next  N statements  define  the  row  of  the  SPLIST 
array  which  correspond  to  this  process.  N being 
the  number  of  pollutants. 

The  next  statement  returns  execution  to  the  PROCESS 
subrout i ne . 

END  statement 
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SUBROUTINE  START  - This  routine  reads  in  the  number  of  sources, 

branches,  parameters,  the  topology  matrix,  and 
selects  the  parameters  to  be  used  from  the  master 
list. 


Variable  Definition 

A (MSORS, 

MBRNC) 


NS 

NB 

NP 


NTOP 

C (MSORS, 
MPARM) 


SPLIST  (MSORS, 
MPARM) 


NPLIST  (MPARM) 


NBRNCH  (MBRNC, 2) 


Statements 


- contains  the  topology  matrix  of  up  to  25 
sources  and  50  branches 

- number  of  sources 

- number  of  branches 

- number  of  parameters 

- number  of  outfalls  from  system 

- contains  parameter  concentration  for  each  of 
25  possible  sources  for  up  to  25  parameters. 
Selected  from  SPLIST 

- master  parameter  concentration  list  passed 
from  the  Process  Model 

- parameter  equivalence  array  indicating  the 
number  of  the  parameter  of  the  master  list 
that  is  associated  with  the  parameters  being 
used 

- contains  the  names  (up  to  8 characters)  of 
up  to  50  possible  branches 


1-  8 subroutine  definition  and  common  area 
9"10  reads  in  number  of  sources,  branches,  parameters,  and 
outfalls  (ns,  NB , NP , NTOP,  respectively) 

11-12  reads  in  the  numbers  of  the  parameters  from  the  master 


A 


I3"I6  reads  in  the  topology  matrix 

17" 1 9 reads  in  the  names  of  the  branches 

20-26  selects  the  appropriate  parameter  concentrations  from 
SPLIST  which  contains  the  information  for  the  master 
list  of  parameters  and  compresses  it  into  array  C 
27-32  initializes  variables  for  COMMON,  calls  CHKl  and  returns 
and  ends 
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I 


L 


SUBROUTINE  CHKl  - This  routine  reads  in  parameter  names,  sample 
frequencies,  number  of  methods  available  for 
analysis,  length  of  survey,  and  outputs  topology 
matrix,  and  source  parameter  information 
Variable  Definition 


A (MSORS,  MBRNC) 

NS 

NP 

NB 

C (MSORS,  MBRNC) 

FLOW  (MSORS) 
LENGTH 

POLN  (MPARM,5) 

I STOP 

XNAME  (MSORS) 

NTEMP  (MPARM) 

NPLIST  (MPARM) 
SAMFRE  (MPARM) 
NMA  (MPARM) 


Statements 


- contains  the  topology  matrix  of  up  to  25  sources 
and  50  branches 

- number  of  sources 

- number  of  parameters 

- number  of  branches 

- contains  parameters  concentrations  for  up  to 
25  sources 

- contains  flow  quantity  from  each  source 

- length  of  the  survey  (days) 

- contains  the  names  of  the  parameters  (up  to 
20  characters) 

- the  number  of  non-competing  parameters 

- contains  the  names  of  each  source  (up  to  ^ 
characters ) 

- contains  the  number  of  methods  available  to 
analyze  each  parameter 

- parameter  equivalence  array 

- contains  sample  frequency  for  each  parameter 

- contains  number  of  methods  available  to  analyze 
each  parameter  (from  master  list) 


1-  9 subroutine  definition  and  COMMON  area 
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10-11  reads  in  length  of  the  survey 

12-13  reads  in  the  number  of  non-competing  parameters 

II4-I7  FORMAT  statements  for  output,  and  page  eject 

18-22  prints  number  of  parameters,  number  of  non-competing 

parameters,  length  of  the  survey,  and  labels 
23"33  reads  in  all  parameter  names,  number  of  methods  and 

sample  frequencies  from  the  mater  list  and  saves  only 
those  that  are  currently  being  used.  Also  the  total 
numbers  of  samples  for  the  survey  are  calculated  based 
on  sample  frequency  and  length  of  the  survey 
3^-37  prints  labels  and  number  of  sources,  branches,  and 
outfalls 

38-41  prints  the  numbers  of  the  sources,  the  names  of  the 
sources,  and  the  rows  of  the  topology  matrix 
42-46  prints  labels  and  the  flow  associated  with  each  source. 

47-49  initializes  counters  and  flag  (for  use  if  more  than  6 
parameters  are  used) 

50  if  more  than  6 parameters  are  present  go  to  30  (statement  60) 
51*57  prints  out  parameter  name,  source  name,  and  the  corre- 
sponding concentrations  for  up  to  6 parameters  at  all 
sources 

58-59  if  more  than  6 parameters  go  to  31  (statement  71),  other- 
wise go  to  the  return  segment. 

6O-6I  there  are  more  than  6 parameters  so  initialize  counters 
the  first  time  this  segment  is  used 
62  if  this  segment  has  been  used  before,  go  to  32  (statement 
65) 
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63*6^  increment  counter  and  transfer  to  33  (statement  5i) 

65  set  I to  the  next  parameter  number 
66-69  determines  the  upper  limit  on  the  parameter  number  (NNP) . 
This  segment  is  used  to  guarantee  that  no  more  than  6 
columns  of  parameters  are  printed  at  a time 
70  transfers  to  print  segment  (statements  51"53) 

71-72  if  all  parameters  have  been  printed,  return,  otherwise, 
go  to  30  (statement  60)  to  select  next  6 parameters 
73"75  CONTINUE,  RETURN,  AND  END  control  statements 
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SUBROUTINE  TOP  - This  subroutine  calculates  the  flows,  concentration  of 
parameters,  mass  per  day  of  conservative  parameters  In 
all  branches  of  the  system.  The  effects  of  treatment 
are  Included.  Also  calls  a subroutine  which  flags 
parameters  to  be  measured  at  specified  measure  points. 
Also  calls  CHK2  for  read  out  of  pertinent  Information. 

Varl?''le  Definitions 


A (MSORS,  MBRNC) 

NS,  NSS 
NP,  NPP 
NB,  NBB 

C (MSORS,  MPARM) 

Y (MBRNC,  MPARM) 

FLOW  (MSORS) 

X (MBRNC,  MPARM) 

EFF  (MBRNC) 

P (MSORS,  MPARM) 

NPLIST  (MPARM) 

PT(25) 


- contains  the  topology  matrix  of  up  to  MSORS  sources 
and  MBRNC  branches 

- number  of  sources  In  a particular  study 

~ number  of  parameters  In  a particular  study 

- number  of  branches  In  a particular  study 

- contains  the  parameter  concentrations  for  NS 
sources  and  NP  parameters 

- output  matrix  of  parameter  concentration  for 
each  of  NB  branches  and  NP  parameters. 

- contains  the  quantity  of  flow  for  each  source 
(MGD) 

- contains  mass  parameter  for  each  of  the  NB 
branches 

- contains  the  flow  for  each  branch 

- contains  the  mass  parameter  for  each  source 
(C  (J,  K)  * FLOW  (J)} 

- lists  the  numbers  (from  the  master  list)  of 
the  parameters  In  this  study 

- branch  parameters  passed  to  treatment 
subroutines 
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EFFT 

MSORS 

MPARM 

MBRNC 

TEST 

CRFLAG 

RANK 

FLGPT 


- branch  flow  passed  to  treatment  subroutines 

- the  program  dimension  for  the  number  of  sources 

- the  program  dimension  for  the  number  of  parameters 

- the  program  dimension  for  the  number  of  branches 

- a flag  that  is  set  if  all  source  flows  are  not 

given 

- a flag  that  is  set  if  TOP  is  called  from  CORRCT 
so  that  read  in  of  treatment  data  is  omitted 

- variable  for  summing  masses  in  FLAG  subroutine 

- flagging  level  for  each  parameter 
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Statements 

1 - 7 

8 


9 - n 

)2 

13-28 


29  - 30 
33  - 35 
36  - 38 

39  - ^*9 


50  - 51 

52  - 53 
5'*  - 56 
57  - 61 

62  - 11*2 


subroutine  definition  and  dimensions 
set  test  flag  = 0 Test  = 1 informs  MAIN  that 
there  is  insufficient  flow  information  to  con- 
t i nue 

zero  out  X(  1 ,J)  array 

flag  to  skip  read  in  if  TOP  is  called  from 
CORRCT 


1 


1 


1 


read  in  treatment  data 
Settle:  AREA 

NEUTRA;  MNEUT,  MENU,  NCAUS,  MACY , 

DOME:  I TYPE,  TEMP,  AREAP,  AREAS,  ABVOL,  TFAREA, 

TFVOL,  MLSS,  YN , R1 , R2 , R3,  R4 , K20,  NN 
read  out  AREA 
change  pH  to 

concentration  C(l,J)  times  FLOW  to  give  P(l,J) 
for  each  source 

Do  loop  to  compute  the  flow  of  each  branch, 

EFF(l).  If  flow  of  any  source  is  not  given, 

print  message  and  return  to  MAIN  with  TEST  = I. 

set  flags  NAC  and  NALK  = 0.  These  will  indicate 

when  alkalinity  and  acidity  are  not  in  NPLIST. 

set  constants  XKl  and  XK2 

determine  numbers  for  NAC  and  NALK 

If  both  NALK  and  NAC  are  in  NPLIST  continue. 

If  not  print  out  message  set  Test  = 1 and  return, 
calculates  the  mass  parameter,  X(l,K),  in  branch  I. 
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63  - 66 
67 


68 


69 


70  - 88 


72 

73  - 75 
76  - 79 


80 

81,87 


88 

89 

90 


91 

92  - 93 


1 

I 

set  variables  and  counters  to  zero 
start  search  through  topology  matrix  to  determine 
branch  inputs 

if  source  J does  not  contribute  to  branch  I,  look 
at  next  source 

if  current  branch  is  not  the  first  branch,  go 
to  180  (statement  89) 

for  the  first  branch  calculate  the  contribution 
to  X(I,K)  for  each  contributing  source  J for  each 
parameter  K.  As  each  source  J is  added,  alkalinity, 
acidity,  and  capacity  factors  CT  are  computed  as 
conservative  parameters. 

Source  {H"*"}  for  computing  CT  is  computed 
cil»  oi2»  computed  for  source  J 

if  source  alkalinity  is  known  CT  and  acidity  are 
computed  from  alkalinity 
for  branch  1 look  at  next  source 
if  acidity  is  given  for  source  J then  CT  and 
alkalinity  are  computed  in  terms  of  acidity, 
for  branch  1 look  at  next  source, 
sets  II  to  one  less  the  current  branch  number 
begins  search  back  through  columns  of  topology 
matrix 

L is  the  current  column  being  checked 
if  source  J appears  in  branch  L go  to  200 
(statement  95)  otherwise  look  at  next  branch  back 
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94 


95  - 101 


96  - 99 


102  - 103 
I 104 

j 105  - 123 

I 

\ 105 

V 

106 

107 

108 

109  - 111 
112 

113  - 117 
118 


if  source  has  not  been  Included  In  some  previous 
branch  go  back  to  130  {statement  70}  and  compute 
the  contribution  to  X (I,K),  including  alkalinity, 
acidity,  and  CT. 

if  source  has  been  included  in  some  previous 
branch  L,  add  the  X (L,KJ  of  that  branch  to  the 
current  X (I,K)  calculation  unless  that  X (L,K) 
has  been  included  in  X (I,K)  when  looking  at  a 
previous  source. 

if  the  amount  from  branch  L has  already  been 
added  to  current  branch,  go  to  230  (statement  104) 
otherwise  increment  NKNT  the  number  of  branches 
that  have  been  added,  and  add  to  the  list 
{NTB  (NKNT)  = L} 

compute  CTl  for  branch  L and  get  composite  CT 
look  at  next  source 

given  the  conservative  values  of  CT  and 
alkalinity  for  a mixture  of  flows  in  branch  I, 

{H'*'}  is  computed  by  iteration  using  flow  weighted 

(H'*'}  as  starting  value. 

approximate  value  of  (IT*’}  for  starting 

molar  value  for  alkalinity 

Composite  CT  with  flow  divided  out 

alkalinity  from  CT  and  approximate  {iT*'} 

test  for  Convergence 

decrement  pH  by  0.1 

recalculate  ALK  and  test  for  convergence 
Increment  pH  by  0.1 


I 

f 


i 

I 

1 

'I 

i 
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1 


119  - 123 
124 

125  - 127 


128  - 131 

132 

133 

134 

135  - 139 
140  - 141 
142 

143  - 148 

149  - 153 
159  - 157 
158  - 159 
160 
I6l 

l62  - 163 


recalculates  ALK  and  test  for  convergence 
this  gives  buffered  pH  for  composite  flow 
of  branch 

checks  to  see  if  the  branch  flow  undergoes  a ' 

new  treatment  in  this  branch.  Note  {the  X(L,K) 
values  previously  added  were  treated  values  if 

1 

appropriate!  if  a treatment  occurs  here  go  to 

280  (statement  128)  else  go  to  next  branch.  l 

) 

prepare  variable  for  treatment  subroutine  | 

i f A ( J , I ) = 2 , then  ca 1 1 NEUTRA  I 

if  A(J, I)  = 3,  then  call  SETTLE  I 

if  A(J,I)  = 4,  then  call  DOME 

calling  statements  i 

■i 

restoring  X(l,K)  after  treatment  ; 

end  of  loop  for  each  branch 

change  mass  variables  back  to  concentrations  j 

(Y  ( I , J) } for  output . 

change  the  {H^}  values  to  pH 

change  A matrix  to  - or  1 entries  for  later  use 

set  NBB  and  NPP 

call  CHKZ  for  read  out 

call  FLAG  subroutine 

RETURN  and  END  statements  .; 
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SUBROUTINE  NEUTRA  - This  subroutine  models  the  neutralization  treatments. 

It  models  acid  or  caustic  neutralization  and  for  each 
allows  one  of  several  neutralizing  agents.  The  fol- 
lowing flags  are  required  and  are  read  in  in  TOP 
subrout i ne 


MNEUT  = 1 

acid  neutral i zat ion 

MNEUT  = 2 

caustic  neutralization 

MENU  = 1 

Ca(0H)2 

MENU  = 2 

CaCO^ 

MENU  = 3 

Na2C0^ 

MENU  = 

NaOH 

NCAUS  = 1 

H^SO^  or  HCI  or  HNO^ 

NCAUS  = 2 

CO2 

MACY  = 1 

H2S0^ 

MACY  = 2 

HCI 

MACY  = 3 

HNO^ 

triable 

Def ini tions 

NPLIST 

(MPARM) 

- lists  the  numbers  from  the  master  list  of  the 
parameters  in  this  study 

PT  (MPARM) 

- branch  parameters  passed  to  treatment  subroutii 

EFFT 

- branch  flow  passed  to  treatment  subroutines 

MPARM 

- the  program  dimension  for  the  number  of  parame 

NP 

- the  number  of  parameters  in  a particular  study 

:atements 

1-  5 

ca 1 1 and 

dimension  statements 

6,  7 

constants 

for  capacity  factor  (CT)  calculation 
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Ik 


8 

PT(I)  to 

9-11 

ci|  . a2  branch  h"^ 

12 

if  caustic  neut ra 1 i ca t i on  go  to  1 80 
for  acid  neutralization  continue 

(statement  42)  , 

13-19 

compute  CT  in  and  ACYIN  with  either 
g i ven 

ALK.  or  ACY. 

in  in 

20-22 

clq  and  aj  for  H =7  after  treatment 

23 

if  acid  neutralization  treatment  is 
NaOH  go  to  104  (statement  25) 

by  Ca(0H)2  or 

24 

if  acid  neutralization  treatment  is 
go  to  105  (statement  29) 

by  CaCO^  or  Na2C0^ 

25 

CT  for  Ca(OH),  or  NaOH 

(104) 

out  2 

26 

ACY  ^ for  Ca(OH).  or  NaOH 
out  2 

27 

if  Ca(0H)2  go  to  106  (statement  33) 

28 

if  NaOH  go  to  107  (statement  79) 

29 

ACY  for  CaCO,  or  Na,C0, 

(105) 

out  3 23 

30 

CT  for  CaCO-  or  Na.CO, 

out  3 23 

31 

if  CaCO^  go  to  108  (statement  88) 

32 

33 

if  Na2C0^  go  to  209  (statement  133) 

(106) 

Ca(OH)  added  to  neutralize  to  pH  = 7 

34 

20%  excess  added 

35-96 

do  loop  to  determine  if  both  Ca  and 

SO,  are  included 
4 

in  study.  If  not  print  message  and 

return  . 

47 

SO,  in  molar  units 
4 

48 

total  Ca  in  molar  units 

49 

active  Ca  in  molar  units 
23 

50-52  ALK,  ACY  out  after  neutralization  (all  multiplied  by 

flow) 

53  calculate  solubility  product  (SOLYP) 

54  if  SOLYP  ^1.32  X lO"^  go  to  20  (statement  69) 

55“57  calculate  CaSO^  precipitate 

58  remaining  SO^ 

59  remaining  Ca 

60-68  determine  dissolved  solids,  total  solids,  suspended  solids, 

and  hardress  return 

69-78  determine  Ca , dissolved  solids,  total  solids,  suspended 

(20)  solids  and  hardness  if  there  is  no  precipitate,  return 

79  Calculate  NaOH  required  for  neutralization 

(107) 

80-87  determine  dissolved  solids,  total  solids,  suspended  solids, 
pH,  ALK,  ACY.  return 

88 

(108)  calculate  CaCO^  required  for  neutralization 

89  201  extra  added 

90-101  do  loop  to  determine  if  both  Ca  and  SO^  are  included  in 

study.  If  not  print  message  and  return. 

102  * EFFT 

(112) 

103-  ALK  and  ACY  after  neutral izat ion  (times  flow) 

104 

105  SOj^  in  molar  units 

106  active  Ca  in  molar  units 

107  solubility  product  SOLYP 

108  SOLYP  f_1.32  X 10‘^  go  to  113  (statement  123) 

109-  determine  CaSO,  precipitate 

111 


24 


112 


} 


1 


» 


113 


114- 

122 


123- 

132 
(113) 

133 
(209) 

134- 

141 


142- 

148 

(180) 

149- 

151 

152 

153- 

157 

158 

159 

160- 

164 


165 
(123) 

166 

167- 

171 


172 

(125) 

173- 

177 


remaining  SO^ 
rema i n i ng  Ca 

determine  dissolved  solids,  total  solids,  suspended  solids 
suspended  solids,  and  hardness.  return 
if  no  precipitate  determine  Ca,  dissolved  solids,  total 
solids,  suspended  solids  and  hardness.  retur 
determine  Na^CO^  for  neutralization 


determine  H^,  ALK,  ACY,  dissolved  solids,  total  solids, 

and  suspended  solids  after  neut ra 1 i zat ion . return 

for  caustic  neutralization  - determine  CT.  and  ALK. 

in  in 

if  either  ALK.  or  ACY.  is  given 
in  in 

after  neutralization,  j,  and  and  for  pH  = 7 

if  CO^  neutra)ization  go  to  127  (statement  178) 

+ 

CT  , ALK  , H ACY  after  neutralization 
out  out 


if  not  H2S0^  neut ra 1 i za t i on  go  to  123  (statement  165) 

H^SO^  required  for  neutralization 

dissolved  solids,  total  solids  and  SO^  after  neut ra 1 zat i on . 
return 

if  not  HCl  neutralization  go  to  125  (statement  172) 


HCl  required  for  neutralization 

dissolved  solids,  total  solids,  and  chlorides  after 

neutralization.  return 

HNOj  required  for  neutralization 

dissolved  solids,  Nitrate/Nitrite,  total  solids  after 
neutralization.  return 


f 
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178  ALKout  for  CO.  neutralization 
(127)  ^ 


179  Clout  for  CO2  neutralization 

180  CO2  required  for  neutralization 

181-  H , ALK,  ACY,  dissolved  solids,  and  suspended  solids 
186 

neut  ra 1 i zat ion 


187  return 


after 


26 


SUBROUTINE  SETTLE  - This  subroutine  models  the  settling  pond  or 


clarifier.  The  clarifier  area  is  read  in  acres 
TOP.  SETTLE  is  called  from  TOP  when  a 3 appears 
in  A(J , I ) . 


Variable  Definitions 
PT(MPARM) 

NPLIST  (MPARM) 

EFFT 

MPARM 

NP 

Statements 


- branch  parameters  passed  to  treatment  subroutine 

- lists  the  numbers  (from  the  master  list)  of  the 
parameters  in  this  study 

- branch  flow  passed  to  the  treatment  subroutine 

- the  program  dimension  for  the  number  of  parameters 

- the  number  of  parameters  in  a particular  study 


1-4 

5-6 

7 

8 
9 

10 
1 1 
12 

13 

14 

15 

16-18 

(15) 


- call  and  dimension  statements 

- compute  the  removal  factor  (RF) 

- starts  DO  loop  for  Settling  solids 
settles  COD  which  is  30^  dissolved 
settles  TOC  which  is  30%  dissolved 

if  KJELDAHL  N is  in  study  go  to  15  (statement  16) 
if  total  solids  is  in  tudy  go  to  25  (statement  22) 
computes  suspended  solids  after  settling 
computes  vol . susp.  solids  after  settling 
computes  turbidity  after  sett'Ing 
go  to  10  for  K 23  (statement  45) 
go  to  F7  rf  ammonia  is  irrcluded  ^statement  20) 


19  go  to  10  if  ammonia  is  not  included 

20  compute  settled  Kjeldahl  N 

(17) 


(statement 
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21 


22-kk 

(25) 


45 

(10) 

46-47 


go  to  10  (next  K)  (statement  45) 

calculates  total  solids  if  any  two  of  the  following  three 
are  given:  total  solids,  dissolved  solids,  suspended  solids 

next  K 

Return  - end 
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SUBROUTINE  DOME 


- This  subroutine  models  the  domestic  or  sanitary 
sewage  treatment.  This  model  includes  a trickling 
filter  with  recycle  or  an  activated  sludge  type 
treatment.  DOME  is  called  from  TOP.  Treatment 
system  data  is  read  in  in  TOP  and  includes: 


TRICKLING  FILTER 
AREAP 
ITYPE 

TFVOL 

TEMP 

R2 

R3 

RA 

TFAREA 

K20 

AREAS 

activated  sludge 

AREAP 

ITYPE 

AREAS 

ABVOL 

MLSS 

YN 

R1 


primary  clarifies  area  in  acre 

0 trick) ing  f i 1 ter 

1 activated  sludge 

trickling  filter  volume  in  acre-feet 
wastewater  temperature  ('^C) 
recycle  from  filter  effluent  (MGD) 
recycle  from  secondary  clarifier  (MGD) 
recycle  to  the  filter  (MGD) 
trickling  filter  area  in  acres 
0.23  for  1”  rock  media 
0.13  for  2 1/2  rock  media 
secondary  clarifies  area  in  acres 

as  above 
as  above 
as  above 

ae rati  an  basin  volume  in  MG 
mixed  liquor  suspended  solids  (Mg/c) 
net  yield  computed  from  plant  log  data 
activated  sludge  recycle  (mGD) 
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Variable  Definitions 


NPLIST  (MPARM)  - lists  the  numbers  (from  the  master  list)  of  the 

parameters  in  this  study 

PT  (MPARM)  - branch  parameters  passed  to  the  treatment 

subrout i ne 

EFFT  - branch  flow  passed  to  treatment  subroutine 

MPARM  - the  program  dimension  for  the  number  of 

parameters 

NP  - number  of  parameters  in  a particular  study 

Statement 

3-6  I n i t i a 1 i ze  da ta 

7-  8 Store  input  water  quality  parameters  in  array  TT 
9 Calculate  removal  factor  for  primary  clarification 

10-12  Calculate  primary  clarifier  effluent  total  solids 

13  Calculate  primary  effluent  TKN 

lA  Calculate  primary  effluent  TSS 

15  Calculate  primary  effluent  VSS 

16  Calculate  primary  effluent  Turbidity 

17  Calculate  primary  effluent  COD 

18  Calculate  primary  effluent  TOC 

19  Calculate  primary  effluent  BOD 

20  Route  to  3 for  trickling  filter  system  or  A for  activated 
sludge  system 

21  Calculate  recirculation  ratio 

22  Calculate  organic  loading 

23  Calculate  recirculation  factor 

2A  Calculate  filter  BOD  removal  factor 


25  Correct  BOD  removal  factor  for  water  temperature 

26  Calculate  filter  effluent  suspended  solids  j 

27  Calculate  filter  effluent  COD 

28  Calculate  filter  effluent  TOC 

29  Calculate  filter  effluent  BOD 

30  Calculate  depth  of  filter  j 

31  Calculate  hydraulic  loading  on  filter 

32  Correct  media  coefficient  for  water  temperature 

33  Calculate  filter  nitrification  factor 

34  Calculate  filter  effluent  nitrates 

35  Calculate  filter  effluent  TKN 

36  Calculate  filter  effluent  ammonia 

37  Calculate  flow  to  secondary  clarifier 

38  Calculate  removal  factor  for  secondary  clarifier 

39  Calculate  secondary  effluent  suspended  solids 

42  Calculate  process  loading  intensity 

43  Calculate  sludge  age 

44  Calculate  flow  to  aeration  basin 

45  Calculate  aeration  basin  effluent  BOD 

46  Calculate  aeration  basin  effluent  COD 

47  Calculate  aeration  basin  effluent  TOC 

49  Calculate  aeration  basin  effluent  suspended  solids 

50  Change  area  of  secondary  from  acres  to  square  meters  | 

51  Calculate  secondary  solids  surface  feed 

52  Calculate  secondary  effluent  suspended  solids 

53  Check  for  nitrification  requirement 

54-57  Initialize  coefficients 


58 


Calculate  aeration  basin  detention  time 

31 


59 


Make  initial  guess  of  Ni trosomonas  concentration 


60-61  Set  Ni trosomonas  limits 

62  Set  iteration  counter  for  Ni trosomonas 

63  Make  initial  guess  of  ammonia  concentration 

6^-65  Set  ammonia  limits 

66  Set  iteration  counter  for  ammonia 

67-70  Make  calculations  on  ammonia  and  Ni trosomonas 

71-82  Test  ammonia  concentrations,  adjust  it  necessary  and 

repeat  steps 

99  Calculate  aeration  basin  effluent  nitrates 

100  Calculate  aeration  basin  effluent  TKN 

101  Calculate  aeration  basin  effluent  ammonia 

103-  Store  changed  parameters  in  array  PT 

10^4 
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SUBROUTINE  FLAG  - This  subroutine  flags  parameters 
at  a source  based  on  fraction  of 
contribution,  and  also  flags  any 
sources  and  branches  selected  by 
Variable  Definition 


for  measurement 
total  mass 
parameters  at 
the  user. 


FLGPT  (MPARM) 


POLN  (MPARM,  5) 

NP 

NS 

NB 

RANK  (MPARM) 

P (MSORS,  MPARM) 
A (MSORS,  MBRNC) 
Y (MBRNC,  MPARM) 

NF 

NX  (50) 

NY  (50) 


NFB 

XNAME  (MSORS) 


- an  array  containing  the  fractional  level  of  total 
mass  for  which  a parameter  will  be  flagged  at  a 
source  (the  first  branch  in  which  a source  appears) 

- an  array  containing  the  names  of  the  parameters 

- number  parameters 

- number  of  sources 

- number  of  branches 

- contains  total  masses  of  parameters 

- contains  mass  of  each  parameter  at  each  source 

- topology  matrix  fup  to  25  sources  and  50  branches) 

- contains  the  concentration  in  each  branch  (50 
possible)  of  each  parameter 

- number  of  flagged  sources  (selected  by  user) 

- contains  the  source  number  (or  branch  number)  for 
" contains  the  parameter  number  for  a selected  flag 

point.  There  is  a one  to  one  correspondence  between 
NX  and  NY  (eg.  NX  (3)  contains  the  source  (or  branch) 
number  for  flagged  point  3,  and  NY  (3)  contains  the 
parameter  number  for  flagged  point  3) 

- number  of  flagged  branches  (selected  by  user) 

- contains  the  name  for  each  source 
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Statements 


]-  8 subroutine  definition,  COMMON  and  DIMENSION  area 

9-10  reads  into  FLGPT  the  minimum  mass  contribution  flagging 

levels  for  the  parameters  (from  the  list  chosen  by  the 
user) 

11-15  prints  out  a title  followed  by  the  name  of  each  parameter 
and  its  flagging  level 

16-18  adds  up  the  mass  for  each  parameter  from  each  source  and 
puts  the  total  in  RANK 

19-28  calculates  the  fraction  of  mass  contributed  by  each 

source  for  each  parameter  (TEMP),  and  if  this  is  greater 
than  or  equal  to  the  flagging  level,  the  topology  matrix 
(a)  is  searched  to  find  the  first  branch  where  the  source 
appears,  and  then  the  corresponding  concentration  in 
the  branch-parameter  concentration  matrix  (Y)  is  flagged 
by  setting  it  negative 

29-31  the  number  of  flagged  sources  (selected  by  the  user)  is 

read  in,  and  if  there  aren't  any,  go  to  80  (statement  51) 

32-3^  each  source  (NX)  and  the  parameter  for  that  source  (NY)  are 
read  i n 

35-^1  a title  is  printed  followed  by  the  name  of  each  source  and 
parameter  flagged 

42-50  searches  the  topology  matrix  to  find  the  first  branch  where 
a source  appears,  then  sets  the  corresponding  entry  in  the 
branch-parameter  concentration  matrix  negative  to  indicate 
the  flag 


34 


51-53 

Sh-Sk 


reads  in  the  number  of  flagged  branches  (selected  by  the 
user)  and  if  this  is  zero,  the  routine  is  ended 
prints  a label  followed  by  the  name  of  each  branch  and 
parameter  flagged,  then  sets  the  appropriate  entry  in 
the  branch-parameter  concentration  matrix  (Y)  negative 
to  indicate  the  flag  for  each  of  the  flagged  points  then 
returns 
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,1 


subroutine  CHK2  - This  subroutine  outputs  parameter  concentration  and 
mass  information 
Variable  Definition 


Y (M6RNC,  MPARM) 


N6 

NP 

NBRNCH  (MBRNC,  2) 


POLIN  (MPARM,  5) 


X (MBRNC,  MPARM) 


EFF  (MBRNC) 


- output  matrix  of  parameter  concentration 
for  each  of  50  possible  branches 

- number  of  branches 

- number  of  parameters 

- name  of  each  branch  (up  to  8 characters) 
for  up  to  50  branches 

- name  of  each  parameter  (up  to  20  characters) 
for  up  to  25  parameters 

- output  matrix  of  parameter  mass  for  each  of 
50  possible  branches 

- flow  in  each  of  50  possible  branches 


Statements 


1-  5 
6 

7-n 

12-15 

16-18 

19 

20-26 


subroutine  definition,  common  area,  and  dimensioning 
sets  flag  to  use  in  determining  which  output  segment  is 
being  done  (concentration  or  mass) 
prints  title  on  page  and  outputs  name  of  branch  and 
flow  in  that  branch 

prints  title  on  page  and  outputs  number  of  branches 
and  parameters 

initializes  counters  and  flag  (for  use  if  more  than  6 
parameters  and  present) 

if  more  than  6 parameters  are  present  go  to  30  (statement  29) 
prints  out  parameter  name,  branch  name,  and  the  corresponding 
concentrations  for  up  to  6 parameters  at  all  branches 
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I 

I 

I 

I 

i 

I 


t 

1 


27  if  more  than  6 parameters  are  present  go  to  3'  (statement  AO) 

28  go  to  5 (statement  A2) 

29*30  there  are  more  than  6 parameters,  so  initialize  counters 
the  first  time  this  segment  is  used 
31  if  this  segment  has  been  used  previously  go  to  32 

(statement  3A) 

32-33  increment  counters  and  transfer  to  33  (statement  20) 

3A  set  I to  the  next  parameter  number 
35-38  determines  the  upper  limit  on  the  parameter  number 

(NNP)  . This  segment  is  used  to  guarantee  that  no  more 
than  6 columns  of  parameters  are  printed  at  a time. 

39  transfers  to  print  segment  (statement  20-26) 

AO-Al  if  all  parameters  have  been  printed  go  to  5 (statement  A2), 

otherwise,  go  to  30  (statement  29)  to  select  next  parameters 
A2-A3  if  both  concentrations  and  masses  have  been  printed  go  to 
60  (statement  52) 

AA-A9  prints  title  and  converts  mass  to  pounds  per  day  (array 
Y is  used  temporarily  for  this  information) 

50-51  sets  the  flag  and  goes  to  55  (statement  16) 

52-57  restores  concentrations  in  array  Y 
58-66  converts  concentration  to  pH 

67-68  RETURN  and  END  control  statements 
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SUBROUTINE  LEVEL  - This  subroutine  provides  the  interface  between  the 


Topological  and  Resource  Models.  It  determines 
from  the  topology  matrix  all  possible  levels  cf 
the  system,  a level  defined  as  a set  of  points 
(branches)  which  completely  characterizes  all 
outfalls  of  the  system.  Once  the  levels  have  been 
found  the  Y matrix  Is  transformed  to  the  P matrix. 


Variable  DEFN 
SUM  (MSORS) 

SUMA  (MSORS) 

SOL  (SI , MBRPl) 


TOP  (100) 

SUMRR 

Y (MBRNC,  MPARM) 

NP 

NS 

NB 

NLEV 

NTOP 

A ( MSORS,  MBRNC) 


- contains  linear  combinations  of  columns  of  the 
topology  matrix 

- contains  previous  contents  of  SUM 

- array  containing  the  branch  number  for  each  level 
(IBN).  sol  (I,  1)  contains  the  number  of  points 
(branches)  in  level  1.  SOL  (1,  2*50)  contains 
the  corresponding  branch  numbers  for  each  point 
in  the  level 

- contains  sequences  of  branches  which  are  being 
treated  as  solution  sequences 

- used  to  sum  the  rows  of  the  topology  matrix 

- output  matrix  of  parameter  concent  rat i on  for 
each  of  50  possible  branches 

- number  of  parameters 

- number  of  sources 

- number  of  branches 

- number  of  levels 

- number  of  outfalls 

- topology  matrix 
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NPLA  (25) 
P (MPARM, 

Statements 
I-  7 
9 

9-15 


16-21 


22 

23 

2A-28 


- number  of  points  in  each  level 
25,  MBRNC)  - parameter  concentration  at  a point  and  a 

level.  P (2,  A,  7)  is  the  concentration  of 
parameter  2 at  point  7 of  level  A. 

subroutine  definition,  COMMON  and  DIMENSION  area 
sets  NP  to  the  number  of  parameters 
checks  to  make  sure  that  the  topology  matrix  has 
all  ones  in  the  last  row  and  if  it  doesn't  prints 
an  error  message  and  executes  a special  return. 

Otherwise,  it  goes  to  35  (statement  16) 
initializes  SUMA,  SUM,  JJ  ( number  of  levels  found), 

I (location  of  the  newly  added  branch  in  a sequence 
in  the  TOP  array),  and  M (the  number  of  the  newly 
added  branch  in  a sequence  stored  in  the  TOP  array) 
sets  the  next  branch  I in  the  TOP  array  (tried 
sequences  array)  to  M 
sets  K to  zero 

tests  the  tried  sequence  contained  in  TOP.  If  the 
new  branch  added  to  the  sequence,  (M  added  to  position 
I of  TOP)  violates  the  solution  criterion  (no  SUM 
entries  can  be  greater  than  1),  control  is  transferred 
to  5 (statement  36).  The  SUM  entries  are  formed  by 
adding  the  column  M of  the  topology  matrix  to  columns 
of  old  column  numbers  contained  in  the  first  1-1 
positions  of  the  TOP  array.  If  any  SUM  (j)  is  0,  K 
is  set  to  1 . 


39 


29 


30-31 

32-35 


36-37 

38-i<0 


41 

42 


43-44 


45 


46 


Tests  the  value  of  K.  If  K equals  zero  a solution 
sequence  has  been  found  and  control  transfers  to  7 
(statement  4l).  If  K is  equal  to  one  the  trial  sequence 
has  not  violated  the  solution  sequence,  but  another 
branch  must  be  added  to  the  sequence, 
stores  the  new  SUM  (j)  totals  on  the  SUMA  array 
adds  one  more  element  on  trial  solution  and  sets  M to 
the  new  branch  to  be  added  equal  to  one  less  than  the 
last  branch.  If  M is  less  than  one,  go  to  47  (statement 
94),  otherwise,  go  to  6 (statement  22)  and  test  the 
new  sequence. 

replace  SUM  with  SUMA  (the  previous  values  of  SUM) 
sets  new  trial  branch  to  be  added  to  trial  sequence 
to  old  trial  branch  minus  1.  If  new  trial  branch  is 
greater  than  or  equal  to  1 control  is  transferred  to 
6 (statement  22)  and  the  new  sequence  is  tested. 
Otherwise,  go  to  47  (statement  94) 
increments  the  number  of  solutions  found  (JJ) 
stores  in  row  JJ  column  1 of  SOL  the  number  of  branches 
in  the  new  solution  sequences 

stores  in  row  JJ  column  2,1+1  the  branch  numbers  of 
the  new  solution  sequences 

if  the  maximum  number  of  solution  sequences  (levels)  has 
been  found  go  to  20  (statement  141) 

if  the  last  branch  in  the  new  solution  sequence  is  not 
equal  to  one,  go  to  5 (statement  36)  and  the  case  is 
handled  as  if  the  last  element  in  the  solution  sequences 

had  been  a violation 
40 


r 


hi  if  only  one  solution  sequence  has  been  found,  go  to  55 
(statement  78) 

A8  set  KK  to  zero 

49-53  determine  the  position  (KK)  of  the  last  common  entry 
in  the  last  two  solution  sequences  found 

54  if  there  are  no  common  entries  go  to  55  (statement  78) 

55  decrements  the  branch  in  the  KK  + 2 position  of  the 
solution  sequence  by  1. 

56  if  IR  is  less  than  or  equal  to  one  go  to  351  (statement  74) 
57-58  set  the  next  branch  in  the  trial  sequence  (m)  to  IR  and 

set  the  number  of  branches  in  the  trial  sequence  to  KK 
59-60  transfer  the  first  KK  elements  in  the  last  solution 
sequence  into  positions  1 to  KK-1  of  TOP 
61  sets  LIHT  to  KK  (number  of  branches  in  trial  sequence) 

62-64  resets  SUMA  and  SUM  arrays  to  zero 

65-70  store  the  sum  of  the  columns  of  topology  matrix 
designated  in  TOP  into  SUMA  and  SUM  arrays 
71  go  to  6 (statement  22) 

72-73  if  KK  (number  of  branches  in  trial  sequence)  is  zero, 

go  to  300  (statement  128)  else  go  to  271  (statement  55) 
74-75  decrement  the  number  of  branches  in  the  trial  sequence 
and  go  to  362  (statement  72) 

76-77  set  the  number  of  branches  in  the  trial  sequence  to 
one  and  go  to  271  (statement  55) 

78-82  determine  KK  when  only  one  solution  sequence  has  been 
found,  and  go  to  300  (statement  128) 
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83  if  only  on  branch  is  in  the  trial  sequence,  go  to  300 
(statement  128) 

84-85  determine  the  last  common  entry  and  if  it  is  less  than 
one,  go  to  651  (statement  92) 

86-87  set  M (the  next  branch  in  the  trial  sequence)  to  IR,  and 
the  number  of  branches  in  the  trial  sequence  to  KK 

88-89  set  TOP  to  the  last  solution  sequence 

90-91  set  LIMT2  to  the  number  of  branches  minus  one  and  go 
to  69  (statement  62) 

92-93  decrement  the  number  of  branches  and  go  to  970 
(statement  83) 

94-97  if  00  solution  sequences  are  found  print  an  error 

message  and  return,  otherwise,  go  to  360  (statement  98) 
98-103  determine  KK  as  the  position  of  the  last  common  entry 
between  the  trial  solution  sequence  whose  last  branch 
is  one  and  the  last  solution  sequence  found. 

104  if  there  is  no  common  entry,  go  to  855  (statement  114) 

105"106  if  the  common  branch  is  less  than  or  equal  to  one,  go 

to  851  (statement  110) 

107-109  save  the  common  branch  number,  the  number  of  branches 
110-111  decrement  the  number  of  branches  and  go  to  the  next 
statement 

112-113  if  there  are  no  branches,  go  to  300  (statement  128), 
otherwise,  go  to  87I  (statement  105) 

114-118  determines  KK  (number  of  branches  in  solution  sequence) 
from  TOP 

119  if  there  is  only  one  branch,  go  to  300  (statement  128) 


120-121  if  the  common  branch  is  one  or  less,  go  to  ?5I  (statement  126) 

122-123  saves  the  common  branch  and  the  location  of  the  newly  added 

branch 

12^-125  sets  L1MT2  to  the  previous  position  and  go  to  69 
(statement  62) 

126-127  decrement  KK  and  go  to  670  (statement  119) 

128-13^  check  to  see  if  any  more  solution  sequences  exist.  If 

not,  go  to  20  (statement  lAl) 

135-136  try  a new  sequence  whose  first  element  (M)  equals  the 

first  branch  in  the  TOP  array  minus  one  and  set  1 to  one. 
137-1^0  reset  SUMA  and  SUM  arrays  to  zero  and  go  to  6 (statement  22) 

lAl-lAS  determines  NMAX  and  NEMAX  as  the  maximum  number  of  branches 

in  a solution  sequence  (level) 

IA6-IA7  defines  DO  loop  counters  L1M2  and  L1M3 

1^*8-155  rearranges  the  rows  of  the  SOL  array  in  increasing 
number  of  branches  in  the  levels 

160  tests  to  see  if  there  is  more  than  one  outfall  from 
the  system.  If  so,  go  to  33  (statement  169) 

161  defines  DO  loop  counter  L1M6  (number  of  solution 
sequences) 

162-168  moves  the  levels  in  SOL  down  and  inserts  the  level  that 

consists  of  only  the  outfall,  and  goes  to  52  (statement  170) 

169  redefines  LIM6 

170-180  calculates  P (I,  J,  K)  array  from  the  Y array  and  level 
information  of  SOL  array 
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l8l-)83  defines  number  of  levels  (NLEV)  , maximum  number  of  branches 


184-185 


in  a level  (NPMAX)  , and  the  number  of  parameters  (NOPL) 
sets  up  NPLA  array 


186-187  RETURN  and  END  control  statements 


SUBROUTINE  RM  - This  subroutine  reads  in  and  prints  out  resource  informa- 
tion, allocates  resources  and  calls  subroutines  when 
needed  to  relieve  violations,  and  prints  out  resource 
allocation  information. 

Va r i ab 1 e Definition 


NL 

NP 

NPLA  (25) 

NPOFEQ 

EQNAME  (MEQ,  5) 
EQTIME  (MEQ) 

NPLIST  (MPARM) 

NMA  (MPARM) 

NTEMP  (MPARM) 

MENAME  (MPARM, 3, 5) 

PCRM  (MPARM,  3) 

PMDATA  (MPARM,  3, 
5) 


number  of  levels 

number  of  parameters 

number  of  points  in  each  level 

number  of  pieces  of  equipment 

names  of  equipment  items 

time  available  for  each  item  of  equipment 

parameter  equivalence  array 

number  of  methods  available  for  each  parameter  (in 
the  current  parameter  list) 

number  of  methods  available  for  each  parameter  in 
the  master  list  that  are  not  in  the  current  list 
the  names  of  the  methods 

the  minimum  acceptable  concentration  for  each  method 
contains  resource  information  for  up  to  3 methods  for 
each  parameter.  The  information  is  stored  as  follows 
first  column  - equipment  item  codes  (up  to  A per 
method),  second  column  - equipment  times  per  sample, 
third  column  - van  space  requirement  for  each  item, 
fourth  column  - analysts'  times,  fifth  column  - cost 
per  sample,  set  up  time,  location  of  analysis,  and 


time  contraint 


VNSP  (MEQ) 

POLN  (MPARM.5) 
CNSTAR  (6) 

IDO  (MPARM,  25) 

PM  (MPARM, 25, 3) 


AMAR  (MPARM,  MP2) 


IVIOL 

USENO  (MEQ) 

SMEQTI  (MEQ) 
SUMM  (6) 

Statements 
1-  10 
1 1 


- van  space  requirement  for  each  item  of  equipment 

- names  of  each  parameter 

- contraints  on  van  space,  analysts'  times,  and  cost 

- if  for  a parameter  I at  level  J all  points  are 
flagged,  then  IDO  (I,  J)  = 1,  otherwise,  it  is  zero 

- array  containing  feasible  method  numbers  (up  to  3) 
for  each  pa.ameter  at  each  level.  If  an  entry  is 
zero,  then  no  method  is  feasible  at  that  /level  as  a 
whole  (although  a method  might  be  feasible  at  a point 
i n a 1 eve  I ) . 

- array  containing  selected  methods.  AMAR  (I,  1) 
contains  the  method  number  for  parameter  I.  AMAR  (I,  2) 
contains  the  level  number  at  which  the  method  is  used, 
and  AMAR  (I,  3"25)  contains  the  number  of  the  method 

at  the  corresponding  point  in  the  level  if  the 
concentration  is  greater  than  zero. 

- violation  code  as  follows:  1 - equipment  violation, 

2 - van  space  violation,  3 • 6 analyst  violation 
(analysts  I - k respectively),  7 ■ cost  violation 

- contains  the  number  of  times  each  item  of  equipment 
is  used 

- sum  of  equipment  time  for  each  item 

- contains  totals  of  van  space,  analysts'  times,  and 
cost 

subroutine  definition,  COMMON  and  DIMENSION  area 
initializes  the  violation  indicator 
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12-  16 


I 


17-  20 
21-  26 

27 

28 

29-  30 

31-  32 
33-  35 
36-  37 

38-  39  , 

40-  73 

74 

75-  86 

87-  89 
90 

91-  92 

93 

94 


prints  out  a lavel  and  the  number  of  points  at  each  level 
reads  in  and  prints  out  the  number  of  pieces  of  equipment 
reads  in  and  prints  out  the  equipment  names  and  times 
available  for  each  item 
initializes  the  counter  for  NPLIST 
DO  loop  to  scan  the  master  list 

determines  the  number  of  methods  available  for  each 
parameter 

DO  loop  to  read  in  all  methods 

DO  loop  to  read  in  resource  information  (PMDATA) 
if  the  parameter  under  consideration  is  in  the  current 
list,  increment  N,  then  look  at  the  next  parameter 
FORMAT  statements  for  reading 

DO  loop  to  look  at  all  paramaters  in  the  current  list 
and  to  print  out  ail  information  read  in  about  them 
call  subroutine  CONCK  to  determine  feasible  methods  and 
make  allocations  for  flagged  points 

read  in  and  print  out  constraints  on  van  space,  analysts' 
times,  and  cost 

initialize  LVIOL  and  skip  output  to  a new  page 
DO  loop  to  allocate  for  all  parameters 

initialize  LREM  (initial  method)  and  LT  (flag  indicating 
method  assignment) 

DO  loop  to  scan  all  levels 

if  all  points  at  the  current  level  are  flagged,  don't 
allocate  any  more  resources  and  look  at  the  next  parameter 


47 


95 

96 

97-  98 


99-102 


103-105 


106-114 


115-118 


119 

120-122 


123 


obtains  the  number  of  points  in  the  current  level 
DO  loop  to  scan  all  methods 

if  there  are  no  more  methods  at  the  current  level  to 
to  the  next  level 

if  a feasible  method  is  found  on  the  first  attempt,  go 
to  2050  (statement  103),  otherwise,  set  flag  to  two 
and  save  the  first  method  and  level  number 
save  the  method  number  in  AMAR,  decrement  the  flag  (LT) 
and  save  the  level  number 

call  the  appropriate  subroutines  to  check  for  violations 
using  the  current  method.  If  there  are  violations  control 
is  transferred  to  2200  (statement  115).  Otherwise,  the 
allocations  are  added  and  allocations  for  the  next  para- 
meter are  started. 

if  this  is  not  the  first  time  through  this  segment,  go 
directly  to  another  method,  otherwise,  save  the  first 
violation  number  to  occur  (and  if  it  is,  an  equipment 
number)  and  then  look  at  another  method 
look  at  another  level 

control  reaches  this  point  if  no  method  at  any  level  was 
found  that  did  not  relieve  a violation.  Resources  are 
allocated  for  the  present  parameter  by  the  first  feasible 
method  and  level 

if  an  equipment  violation  could  not  be  relieved,  and 
the  violation  is  on  the  same  equipment  item  as  the  initial 
violation,  go  to  4000  (statement  132) 


124 


if  a violation  could  not  be  relieved  and  it's  the  same 


125-126 

127-128 

129-130 

131 

132-136 

1 37-140 
141-145 
146-149 

150-151 

152 

153 


violation  as  the  initial  violation,  go  to  the  statements 
dependent  on  the  type  of  violation 

attempts  to  relieve  violation  on  equipment.  If  successful 
go  to  1000  (statement  149).  If  not  go  to  4000  (statement  132) 
attempts  to  relieve  violation  on  van  space.  If  successful 
go  to  1000  (statement  149).  If  not  go  to  4100  (statement 

137). 


attempts  to  relieve  violation  on  analyst  time.  If 
successful  go  to  1000  (statement  149).  If  not  go  to 
4200  (statement  1 41 ) . 

attempts  to  relieve  violation  on  cost.  If  successful 
go  to  1000  (statement  149).  If  not  go  to  4300  (statement 
146). 

print  out  violation  information  and  save  violation  number 
and  equipment  number,  and  go  to  1000  (statement  149) 
print  out  violation  information  and  save  violation  number, 
and  go  to  1000  (statement  149). 

print  out  violation  information  and  save  violation 

number,  and  go  to  1000  (statement  l49) 

print  out  violation  information  and  save  violation 

number,  and  look  at  the  next  parameter 

calls  subroutines  to  print  out  measure  information  and 

to  allocate  time  for  sampline  and  flow  measurement 

initializes  flag  point  counter 

DO  loop  to  look  at  all  parameters 
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154-159 

160-162 

163 

164 

165-170 

171-175 

176-182 

183-198 

199-205 

206-207 


calls  subroutine  to  print  out  parameter  allocation 
informat  ion 

prints  out  a label  and  calls  a subroutine  to  obtain 
equipment  usage  information 

DO  loop  to  look  at  all  possible  equipment  items 
if  an  item  of  equipment  has  not  been  used,  go  to  5200 
(statement  175) 

initializes  variables  and  determines  if  an  equipment 
time  violation  occurred,  and  the  amount  of  time  an 
i tern  was  used 

calls  subroutine  to  convert  minutes  to  hours  and  minutes, 
and  prints  out  equipment  usage  information,  then  looks 
at  the  next  equipment  item 

prints  out  van  space  usage  information  (van  space  used, 
if  a violation  occurred,  and  the  amount  of  the  violation) 
prints  out  the  grand  total  of  each  analysts'  time,  of  a 
violation  occurred,  and  the  amount  of  the  violation 
prints  out  the  grand  total  of  cost,  if  a violation 
occurred,  and  the  amount  of  the  violation 
RETURN  and  END  control  statements 


SUBROUTINE  CONCK  - This  subroutine  determines  feasible  methods  at  all 


levels  for  all  parameters,  and  allocates  resources  to 
flagged  points. 

Variable  Definition 


NP 

POLN  (MPARM,  5) 
MENAME  (MPARM,  3,  5) 
FP  (i*00,  4) 


IFPT 

IDO  (MPARM,  25) 


NL 

NMA  (MPARM) 

NPLA  (25) 

PC  (MPARM,  25,  MPARM) 


PChM  (MPA“M,  3) 

PM  (MPARM,  25,  3) 


- number  of  parameters 

- names  of  parameters 

- names  of  methods  (up  to  3 per  parameter) 

- flagged  point  information,  FP  (l,  I)  contains 
the  parameter  number  for  flagged  point  I. 

F’^  (l,  2)  contains  level  number,  FP  (l,  3) 
contains  point  number,  and  FP  (l,  4)  contains 
the  method  number 

- number  of  f ta-gqg<  points 

- if  for  a parameter  I , at  a level  J,  all 
points  are  flagged,  IDO  (I,  J)  * 1,  otherwise 
it  is  zero. 

- number  of  levels 

- contains  number  of  methods  available  to  analyze 

- contains  the  number  of  points  at  each  level 

- contains  the  parameter  concentration  information. 
PC  (I,  J,  K)  contains  the  concentration  of 
parameter  I at  point  K of  level 

- contains  the  minimum  acceptable  concentration 
for  each  method  (3  possible)  for  each  parameter 

- array  containing  feasible  method  numbers  (up 
to  3)  for  each  parameter  at  each  level 
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IBN  (51  , MBRPl ) 


NALOW  (MPARM,  MBRNC) 


SMEQTI  (MEQ) 

SAMFRE  (MPARM) 

EQUSED  (MPARM, .MEQ) 

EU  (MEQ) 

SUMM  (6) 

LENGTH 

NSET  (MPARM,  3) 

PMDATA  (MPARM,  3,  5) 


array  containing  the  branch  numbers  for  each 
level.  IBN  (I,  1)  contains  the  number  of 
points  (branches)  in  level  I.  IBN  (I,  2-50) 
contains  the  corresponding  branch  numbers  for 
each  point  in  1 eve  I I . 

contains  flags  to  determine  if  a parameter 
is  measured  at  a branch.  If  NALOW  (I,  J) 
is  a one  then  parameter  I is  measured  at 
branch  J,  zero  if  it  is  not  measured, 
sum  of  equipment  time  for  each  item 
the  sample  frequency  of  each  parameter 
number  of  times  a piece  of  equipment  has 
been  used  for  a particular  parameter 
total  number  of  times  each  item  of  equipment 
is  used 

contains  sums  of  van  space,  analysts'  times, 
and  cost  for  previous  allocations 
length  of  the  survey 

array  to  indicate  if  set  up  time  for  a given 
method  has  already  been  added  to  the  total 
array  containing  resource  information  for  up 
to  3 methods  for  each  parameter.  The  informa- 
tion is  stored  as  follows:  first  column  - equipment 

codes  (up  to  I4  items  for  each  method),  second 
column  - equipment  time  per  sample,  third 
column  - van  space  for  each  item,  fourth  column- 
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1 


analysts'  times,  fifth  column  - cost  per  sample, 
set  up  time,  location  of  analysis,  and  time 
constra  i nt . 


Statements 
1-  6 

7 

8 

9 
10 
1 1 


12 

13-U 

15 

16 


17 

18-20 

21 

22 

23 

2U 

25 


subroutine  definition,  COMMON,  and  DIMENSION  area 
DO  loop  to  look  at  all  parameters 

sets  counter  to  zero  (MNCK  indicates  if  a feasible  method 
has  been  found) 

DO  loop  to  look  at  all  levels 
sets  method  counter  (MC)  to  zero 

determines  the  number  of  methods  available  for  the  current 
parameter 

DO  loop  to  look  at  all  methods 

DO  loop  to  check  all  points  at  the  current  level 
if  the  concentration  is  zero,  go  to  AOO  (statement  17) 
if  the  concentration  is  less  than  the  minimum  acceptable 
for  the  current  method  go  to  300  (statement  21) 
look  at  the  next  point  in  the  level 

increment  MNCK  to  indicate  a feasible  method  has  been 
found,  increment  the  method  counter  and  store  the  number 
of  feasible  method 
look  at  the  next  method 
look  at  the  next  level 

if  no  feasible  method  was  found,  allocate  the  first 

method  anyway 

look  at  the  next  parameter 

initialize  the  flag  point  counter  to  zero 


J 
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26-27  DO  loops  to  look  at  all  parameters  and  all  levels 

28-30  initialize  flagged  point  indicator  array,  flagged  points 

at  a level  and  feasible  method  indicator 
31-32  DO  loop  to  look  at  all  points  in  the  current  level 

33  if  a point  is  not  flagged  go  to  700  (statement  72) 

3^  if  a method  is  available  go  to  750  (statement  kk) 

35-36  determine  the  number  of  methods  available  and  look  at 

w 1 1 of  them 

37“38  if  the  method  is  feasible  go  to  735  (statement  ^2), 

otherwise,  look  at  the  next  method 
39  set  LN  to  indicate  no  feasible  method  was  found 

^0-^1  since  no  method  was  feasible  use  the  first  method  anyway 

and  go  to  750  (statement  kk) 

42-43  save  the  feasible  method  number  and  set  LN  to  indicate 

a feasible  method  was  found  at  a point 
44-46  determine  the  branch  that  the  levei  and  point  correspond 

to,  and  if  a measurement  has  already  been  made  there,  go 
to  1000  (statement  67),  otherwise,  set  NALOW  to  indicate 
that  a measurement  will  have  been  made 

47  obtain  the  method  number 

48  DO  loop  to  add  all  equipment  times,  van  space,  analysts' 
times  and  cost 

49-53  if  no  equipment  item  is  used,  go  to  910  (statement  56), 

otherwise,  add  new  equipment  time  to  previous  total,  and 
increment  the  number  of  times  the  item  is  used 
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J 


5^-55  if  the  item  has  been  used  previously,  go  to  910  (statement 

56),  otherwise,  add  van  space  requirement  to  previous 
total 

56-58  if  set  up  time  has  already  been  added,  go  to  900  (statement 

59),  otherwise,  add  set  up  time  to  the  appropriate 
analysts'  total,  then  set  the  array  to  indicate  that 
set  up  time  has  now  been  added 

59-60  add  sample  time  requirement  to  the  analysts'  total  and 

close  the  loop  (statement  ^8) 

61  add  new  cost  to  the  previous  total 

62-66  increment  the  flagged  point  counter  and  save  the  parameter 

number,  level  number,  point  number,  and  method  number 

67  if  no  feasible  method  existed  go  to  925  (statement  70) 

68-69  if  a feasible  method  existed  only  at  a point,  reset  PM 

to  zero  and  go  to  9^0  (statement  71) 

70  if  the  level  is  not  the  first  level,  reset  PM  to  zero 

71-72  increment  the  number  of  flagged  points  at  a level  and 

look  at  the  next  point 

73  if  all  points  at  a level  are  flagged,  set  IDO  to  one 

74  look  at  the  next  level 

75  look  at  the  next  parameter 

76-77  RETURN  and  END  control  statements 
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SUBROUTINE  TO  - This  subroutine  is  used  to  zero  out  temporary  allocations 
for  equipment  times. 

Variable  Definition 

TEMP  {k , 2)  - contains  equipment  number  and  total  time  allocated 

for  that  item  for  the  method  under  consideration 
(up  to  k items  per  method) 

Statements 

1-2  subroutine  definition  and  dimension  statement 

3-6  zeros  out  TEMP 

7-8  return  and  end 
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SUBROUTINE  EOCHEC  - This  subroutine  determines  if  any  equipment  time 

violations  occur  for  a method  under  consideration. 


Variable  Defini 
I 

J 

L 

MC 

NPLAJ 
I VIOL 
lEQN 

PC  (MPARM,  25 
MBRNC) 

AMAR  (MPARM, 
MP2) 


PMDATA  (MPARM, 
3,'4,5) 


TEMP  (1*,2) 


t ion 

- current  parameter  number 

- current  level  number 

- current  method  number 

- not  used 

- violation  code 

- violation  code 

- equipment  number 

, -parameter  concentration  at  a level  and  a point  (eg.  PC 
(2,  A,  7)  is  concentration  of  parameter  2 at  point  7 of 
level  A) 

- array  containing  feasible  methods.  AMAR  (I,  1)  contains 
the  method  for  parameter  I.  AMAR  (I,  2)  contains  the 
level  that  the  method  is  used  at,  and  AMAR  (I,  3"25) 
contains  the  number  of  the  method  if  it  is  used  at 

the  corresponding  branch  in  the  level 

- array  containing  resource  information  for  up  to  3 methods 
for  each  parameter.  The  information  is  stored  as  follows: 
first  column  - equipment  codes  (up  to  A items  for  each 
method),  second  column  - equipment  time  per  sample, 

third  column  - van  space  for  each  item,  fourth  column  - 
analysts'  time,  fifth  column  - cost  per  sample,  set  up 
time,  location  of  analysts'  and  time  constraint 

- contains  equipment  number  and  total  time  allocated 
for  that  item  for  the  method  under  cons i oerat i on 
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SAMFRE  (MPARM)  - sampling  frequency  for  each  parameter  (for  the  total 


SMEQTI  (MEQ) 
EQTIME  (MEQ) 
Statements 
1-  A 
5 

6-  9 


10 

11-22 


survey)  j 

- sum  of  equipment  time  for  each  item  ] 

- total  time  available  for  each  item  of  equipment  j 

subroutine  definition,  COMMON  and  DIMENSION  area 
DO  loop  definition,  to  check  all  points  at  the  level 
under  consideration 

if  parameter  concencrat ion  at  point  K of  level  J is 
greater  than  zero,  go  to  150  (statement  11),  otherwise 
set  the  corresponding  branch  entry  in  AMAR  to  zero  and 
check  the  next  branch 

sets  the  corresponding  branch  entry  in  AMAR  to  the  method 
under  consideration 

selects  the  number  of  the  equipment  item  (if  no  equipment, 

I EQN  = 0 so  skip  the  time  check),  calculates  the  time 
used  for  the  method  under  consideration,  and  checks  if  this 
will  violate  a time  constraint.  If  it  does  then  the 
error  flag  is  set  (IVIOL  ” 1)  and  the  special  return  is 
used.  Otherwise,  a normal  return  is  executed.  This  is 
done  for  each  of  four  possible  items  of  equipment. 
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SUBROUTINE  VSCHEC  - This  subroutine  checks  for  van  space-cons t ra i nt 

V i o 1 a t ions . 

Var i able  Def i n i t i on 


I 

L 

IVIOL 

SUM 

PMDATA  (MPARM,  3,  5)- 


EU  (MEQ) 
SUMM  (6) 

CNSTAR  (6) 


current  parameter  number 
current  method  number 
violation  code 

van  space  used  in  current  method 
array  containing  resource  information  for  up 
to  3 methods  for  each  parameter.  The  informa- 
tion is  stored  as  follows:  first  column  - equipment 

item  codes,  second  column  - equipment  time  per 
sample,  third  column  - van  space  for  each 
item,  fourth  column  - analysts'  times  per 
sample,  fifth  column  - cost  per  sample, 
set  up  time,  location  of  analysis,  time  constraint 
number  of  times  each  item  of  equipment  is  used 
contains  sums  of  van  space,  analysts'  times, 
and  cost  for  previous  allocations 
contains  constraints  for  van  space,  analysts' 
times,  and  cost 


Statements 


1-  3 subroutine  definition,  COMMON  and  DIMENSION  area 

4 sets  current  van  space  allocation  to  zero 

5-11  computes  van  space  needed  for  current  method.  If  an  item 

of  equipment  has  been  used  previously  (EU  (lEQN)  greater 
than  1),  it  is  not  added  to  the  current  van  space 
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12-16 


determines  if  a van  space  violation  occurs.  If  not  a 


normal  return  is  executed.  Otherwise,  the  violation 
code  is  set  and  a special  return  is  executed. 
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SUBROUTINE  CNCHEC  - This  subroutine  determines  if  any  analyst  time 

violations  or  cost  violations  occur. 

Variable  Definition 

ICONNO  - analyst  number  (or  row  number  of  cost  information) 

I - current  parameter  number 

J - current  level  number 

L - current  method  number 

ICP  - column  number  of  analyst  (or  cost)  information 

IVIOL  - violation  code 

TOT  - total  of  an  analysts'  time  (or  cost)  for  current 

method 

NPLA  (25)  ■ number  of  points  in  each  level 

NSET  (MPARM,  3)  ■ array  to  indicate  if  set  up  time  for  a given 

method  has  already  been  added  to  the  total 
LENGTH  - length  of  the  survey 

AMAR  (MPARM,  MP2)  - array  containing  possible  methods.  AMAR  (l,  1) 

contains  the  method  for  parameter  I.  AMAR  (I,  2) 
contains  the  level  that  the  method  is  used, 
and  AMAR  (I,  3~25)  contains  the  number  of  the 
method  at  the  corresponding  point  in  the  level 
if  the  concentration  is  greater  than  zero. 

SAMFRE  (MPARM)  - sampling  frequency  for  each  parameter 

SUMM  (6)  - contains  sums  of  van  space,  analysts'  times 

and  cost  for  previous  allocations 

CNSTAR  (6)  - contains  constraints  of  vin  space,  analysts' 

times  and  cost 
61 
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PMDATA  (MPARM,  3.^.5)  " array  containing  resource  information  for 


r 


up  to  3 methods  for  each  parameter  (see 
variable  definition  for  SUBROUTINE  EQCHEC) 

Sta  tements 

1-  h subroutine  definition,  COMMON  and  DIMENSION  area 

5 sets  current  time  total  (or  cost)  to  zero 

6 sets  NPLAJ  to  the  number  of  branches  in  level  J 

7 if  cost  is  being  calculated  go  to  20  (statement  14) 

8 if  set  up  time  has  already  been  added  for  this  method 
go  to  20  (statement  14) 

9 if  analyst  ICONNO  is  used  for  this  method,  calculate  the 
set  up  time 

10-13  for  each  branch  of  the  level  for  which  the  current  method 

is  used  calculate  analysts'  time  (or  cost)  and  add  it  to 
the  total 

14  if  analysts'  times  are  being  considered  go  to  15  (statement  18) 

15  if  no  cost  was  added,  return  without  checking  for  a violation 

16-17  if  no  cost  violation  return  normally,  otherwise,  go  to 

50  (statement  21) 

18  if  no  time  was  added,  return 

19-20  if  no  time  violation  return  normally,  otherwise,  go  to 

50  (statement  21) 

21-24  if  a violation  has  occurred,  set  the  proper  violation 

code  and  execute  a special  return  and  END  control  statement 
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SUBROUTINE  ADD  - This  subroutine  adds  the  current  allocations  of  equipment 
time,  van  space,  analysts'  times,  and  cost  to  the  previous 


tota 1 s . 

Variable  Def i n i t ion 
I 

J 

L 

NPLA  (25) 

PC  (MPARM,  25, 

MBRNC) 


PMDATA  (MPARM,  3, '*,5) 


SMEQTI  (MEQ) 

SAMFRE  (MPARM) 

EQUSED  (MPARM,  MEQ) 

EU  (MEd) 


current  parameter  number 
current  level  number 
current  method  number 
number  of  branches  in  each  level 
parameter  concentration  at  a level  and  a point 
(eg.  PC  (2,  4,  7)  is  the  concentration  of  para- 
meter 2 at  a point  7 of  level  4) 
array  containing  resource  information  for  up 
to  3 methods  for  each  parameter.  The  informa- 
tion is  stored  as  follows;  first  column  - 
equipment  codes  (up  to  4 items  for  each  method), 
second  column  - equipment  time  per  sample, 
third  column  - van  space  for  each  item,  fourth 
column  - analysts'  time,  fifth  column  - cost 
per  sample,  set  up  time,  location  of  analysis 
and  time  constraint 
sum  of  equipment  time  for  each  item 
sampling  frequency  for  each  parameter  (for  the 
total  survey) 

number  of  times  a piece  of  equipment  has  been 
used  for  a particular  parameter 
total  number  of  times  each  item  of  equipment 
is  used 


3 

i 


I 


i 
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I 


f 


SUMM  (6) 


NSET  (MPARM,  3) 


Statements 


- contains  sums  of  van  space,  analysts'  times, 
and  cost  for  previous  allocations 

- array  to  indicate  if  set  up  time  for  a given 
method  has  already  been  added  to  the  total 


7 

8 
3 

10 


11-12 


13 

lA-15 

16 


17 

18 

19-20 


21 

22 

23-25 


subroutine  definition,  COMMON  and  DIMENSION  area 

sets  NPLAJ  to  the  number  of  branches  in  level  J 

00  loop  to  look  at  each  branch  in  level  J 

if  parameter  has  zero  concentration  in  a given  branch  do 

not  add  anything  to  the  totals;  check  the  next  branch 

DO  loop  to  look  at  all  possible  equipment  items  used 

and  analysts'  used 

if  no  equipment  is  used  don't  add  any  equipment  usage 
transfer  to  300  (statement  I8) 

add  the  time  ^ the  appropriate  equipment  time  total 
i ncrement  the  number  of  times  a piece  of  equipment  is  used 
if  an  item  of  equipment  has  been  used  before,  go  to  300 
statement  (statement  l8)  to  avoid  add i ng  van  space  again 
add  van  space  requirement  for  equipment  item 
add  analysts'  time  ^ total 

if  set  up  time  has  been  allocated  previously  go  to  200 
(statement  21)  to  avoid  add i ng  it,  otherwise,  if  an  analyst 
is  used  for  this  method,  then  add  the  set  up  time 
end  of  DO  loop 
adds  cost  ^ total 

end  of  DO  loop  and  RETURN  and  END  control  statements 
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SUBROUTINE  SET  - This  subroutine  sets  up  the  array  of  methods  used  at 
a branch  in  a level  for  a given  parameter. 

Variable  De i i n i t ion 

I - current  parameter  number 

J - current  level  number 

L - current  method  number 

PC  (MPARM,  25,  - parameter  concentration  at  a point  and  a level 

MBRNC) 

(eg.  PC  (2,  k,  7)  is  the  concentration  of  a 
parameter  2 at  point  7 of  level 

A (MPARM,  MP2)  - array  containing  possible  methods  (AMAR) . 

A (I,  1)  contains  the  method  for  parameter  1. 

A (1,  2)  contains  the  level  number  that  the 
method  is  used,  and  A (l,  3*25)  contains  the 
number  of  the  method  at  the  corresponding 
point  in  the  level  if  the  concentration  is 
greater  than  zero. 

N (25)  - number  of  points  in  each  level  (NPLA) 

Statements 

I-  3 subroutine  definition  and  DIMENSION  area 

k sets  A (1,  1)  to  the  current  method  number 

5 sets  A (1,  2)  to  the  current  level  number 

6 sets  NP  to  the  number  of  points  in  level  J 

7 DO  loop  to  check  every  point  in  level  J 

8-10  if  the  concentration  at  point  K in  level  J is  greater 

than  zero  go  to  kO  (statement  11),  otherwise,  set  A 
(1,  K + 2)  to  zero  and  check  the  next  point. 
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1 


11-12 

13-'^ 


sets  the  corresponding  point  in  array  A to  the  current 
method  number  and  check  the  next  point 


RETURN  and  END  control  statements 
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SUBROUTINE  PIEQCH  - This  subroutine  searches  past  assignments  to  attempt 

to  alleviate  constraint  violations 

Variable  Definition 


I STOP 
lEQNR 

NOV 

AMAR  (MPARM,  MP2) 


PM  (MPARM,  25,  3) 

NPLA  (25) 

PMDATA  (MPARM,  3,  , 5)- 


NL 


current  parameter  number 

number  of  noncompeting  parameters 

code  of  equipment  items  for  which  a violation 

has  occurred 

violation  code 

array  containing  allocated  methods.  AMAR  (l,  I) 
contains  the  method  number  for  parameter  I. 

AMAR  (1,  2)  contains  the  level  number  at  which 
the  method  is  used,  and  AMAR  (l,  3*25)  contains 
the  number  of  the  method  at  the  corresponding 
point  in  the  level  if  the  concentration  is 
greater  than  zero. 

array  containing  feasible  method  numbers 
(up  to  3)  for  each  parameter  and  each  level 
number  of  points  at  each  level 
array  containing  resource  information  up  to 
3 methods  per  parameter.  The  information  is 
stored  as  follows:  first  column  - equipment 

codes  (up  to  k items  for  each  method),  second 
column  - equipment  time  per  sample,  third  column  - 
van  space  for  each  item,  fourth  column  - analysts' 
times,  fifth  column  - cost  per  sample,  set  up 
time,  location  of  analysis,  and  time  constraint 
number  of  levels 
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Statements 


I-  5 subroutine  definition,  COMMON  and  DIMENSION  area 

6-  7 prints  out  violation  information 

8 sets  IMO  to  the  current  parameter  number 

9"10  looks  at  the  next  lower  parameter  and  if  it  is  a 

noncompeting  parameter,  RETURN  1. 

11-12  obtains  the  method  number  (LR)  and  the  level  number  (jR) 

13-1^  saves  the  method  and  level  numbers  first  tried 

15"16  prints  out  method  information 

17  if  the  violation  is  on  analysts'  time,  go  to  200  (statement  24) 

18  if  the  violation  is  on  van  space  or  cost,  go  to  300 

(statement  25),  otherwise,  the  violation  was  on  equip'^ent 
t ime 

19-23  searches  to  determine  if  the  current  method  used  the 

equipment  item  on  which  the  violation  occurred  and  if  it 
does,  go  to  300  (statement  25).  If  it  doesn't  go  to  50 
(statement  9) . 

24  if  the  current  method  does  not  use  the  analyst  for  which 
there  was  a violation,  go  to  50  (statement  9)- 

I 

25  call  subroutine  SUBT  j 

26-28  determines  if  the  current  method  is  a feasible  method  and  | 

I 

if  it  is,  go  to  500  (statement  29)  | 

• i 

29-31  if  the  method  is  the  last  (third)  feasible  method,  go  to  I 

550  (statement  32),  otherwise,  set  MC  to  zero  and  increment 
the  level  number 

32  increment  MC 

I 

I 


J 
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33  DO  loop  to  search  all  rernaining  levels 

sets  NPLAJ  to  the  number  of  points  in  level  J 
35  DO  loop  to  look  at  all  feasible  methods 

36-37  if  there  is  not  a feasible  method  for  the  current  level  go 

to  690  (statement  52) 

38-51  call  subroutines  to  see  if  new  method  relieves  violation, 
and  if  it  does  call  ADD,  set  the  violation  code  to  zero, 
and  use  RETURN  2.  Otherwise,  check  another  method. 

52-53  set  MC  to  one  and  check  another  level 

54-55  if  violation  was  not  removed,  add  in  time  for  previously 

allocated  method,  and  check  another  parameter 

56-57  RETURN  and  END  control  statements 


I 

I 

I 

f 

I 

I 

i 

I 
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SUBROUTINE  SUBT  - This  subroutine  subtracts  the  current  allocations 
of  equipment  time,  van  space,  analysts’  times,  and 
cost  from  the  previous  totals 
Variable  Definition 


same  as  SUBROUTINE  ADD 
Statements 

same  as  SUBROUTINE  ADD,  except  replace  add  w i th  "subtract", 
to  with  "from",  i ncrement  w i th  "decrement" 
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SUBROUTINE  INFORM  - This  subroutine  prints  out  the  name  of  every 


branch  and  the  names  of  all  the  parameters 
measured  at  each  branch. 

Var i able  Definition 


AMAR  (mPARM,  MP2) 


I BN  (51  , MBRPI) 


NPLA  (25) 

POLN  (MPARM,  5) 
FP  (AOO,  k) 


IFPT 

NBRNCH  (MBRNC,  2) 
BRN  (MBRNC) 


- array  containing  selected  methods.  AMAR  (l,  1) 
contains  the  method  for  parameter  I.  AMAR  (I,  2) 
contains  the  level  number  at  which  the  method  is 
used,  and  AMAR  (l,  3"25)  contains  the  number  of 
the  method  at  the  corresponding  point  in  the 
level  if  the  concentration  is  greater  than  zero 

- array  containing  the  branch  numbers  for  each 
level.  I BN  (l,  1)  contains  the  number  of  points 
(branches)  in  level  I.  I BN  (l,  2-50)  contains 
the  corresponding  branch  numbers  for  each  point 
in  the  I eve  1 . 

- contains  the  number  of  points  (branches)  in  each 
1 eve  1 

- contains  the  name  of  each  parameter 

- contains  flagged  point  information.  FP  (1,1) 
contains  the  parameter  number  for  flagged  point  I. 
FP  (l,  2)  contains  the  level  number  for  flagged 
point  I,  FP  (l,  3)  contains  the  point  number 

for  the  level  indicated,  and  FP  (l,  k)  contains 
the  method  used. 

- the  number  of  flagged  points 

- contains  the  name  of  each  branch 

- contains  the  number  of  parameters  measured  at 


each  branch 

71 


number  of  parameters 
number  of  branches 


Statements 


17-18 


28-30 


subroutine  definition,  COMMON  and  DIMENSION  area 
zeros  out  the  array  BRANCH  and  BRN 
prints  a label 

DO  loop  to  look  at  all  parameters 

sets  NPTS  to  the  number  of  points  in  the  level  obtained  from  AMAR 
DO  loop  to  look  at  all  points 

if  a parameter  is  not  measured  at  a point,  go  to  100  (statement  2 1 
sets  N to  the  branch  number  and  increments  the  number  of  parameters 
measured  at  that  branch 

places  the  parameter  number  in  BRANCH  and  closes  the  DO  loop 
DO  loop  to  look  at  all  branches 

Initializes  KOUNT  (a  flag  used  to  determine  if  any  measurements 
were  made  at  a branch) 

DO  loop  to  look  at  all  parameters 

prints  out  the  branch  name  the  first  time  through  the  loop 
if  no  more  parameters  were  allocated  at  a branch,  go  check 
the  flagged  points 

set  KOUNT  to  indicate  a measurement  was  made  at  a branch 
prints  the  name  of  the  parameter  measured  and  continues  the 


DO  1 oop 

31  if  there  are  no  flagged  points,  go  to  211  (statement  kO) 

32  DO  loop  to  look  at  all  flagged  points 

33-35  obtains  parameter  number  (NSUB) , level  number  (MM),  and  point 
number  (NN) 


36-35  obtains  the  branch  number,  and  if  it  is  the  same  as  the 

branch  under  consideration,  print  the  name  of  the 
parameter,  and  set  the  flag  (KOUNT)  to  indicate  that  a 
measurement  was  made  at  that  branch,  then  look  at  the 
next  flagged  point 

^0-^1  if  no  measurements  were  made  at  a branch,  print  this 

information  and  look  at  the  next  branch 
^2-50  CONTINUE  statements,  FORMAT  statements,  and  RETURN,  and 

END  control  statements 
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SUBROUTINE  SAMPLE  - This  subroutine  allocates  analyst  time  (c I ass i f i cat  ion  I) 


for  taking  samples  and  flow  measurements.  The  method 
used  for  sampling  at  each  branch  is  determined  by  the 
user,  and  allocation  of  time  is  based  on  the  highest 
number  of  samples  required  at  a branch.  Flow  meas- 
urement is  also  determined  by  the  user. 


Variable  Def i n i t i on 
SAMFRE  (MPARM) 
SUMM(6) 


IFPT 

FP  (400,  4) 


I BN  (51  , MVRPl) 


NBRNCH  (M8RNC,2) 

NP 

NB 

LENGTH 


- sample  frequency  for  each  parameter 

- contains  allocations  of  van  space,  analysts'  times 
and  cost.  SUMM  (2)  contains  time  allocation  for 
analyst  c 1 ass i f i cat i on  1. 

- the  number  of  flagged  points 

- contains  flagged  point  information.  FP  (l,  1) 
contains  the  parameter  number  for  flagged  point  1. 

FP  (I,  2)  contains  the  level  number,  FP  (l,  3)  contains 
the  point  number,  and  FP  (l,  4)  contains  the  method 
number . 

- array  containing  the  branch  numbers  for  each  level. 

I BN  (I,  1)  contains  the  number  of  points  (branches) 
in  level  I.  IBN  (I,  2-50)  contains  the  corresponding 
branch  numbers  for  each  point  in  the  level. 

- contains  the  name  for  each  branch 

- number  of  parameters 

- number  of  branches 

- length  of  the  survey 
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BRANCH  (MBRNC,  MPARM) 

NSAMPL  (10,5) 

SAMPTI  (10) 

SETUP  (10) 

NROUT  (MBRNC) 

NFLOW  (MBRNC) 


- contains  the  number  of  the  parameters 
measured  at  each  branch 

- contains  the  names  of  the  sampling  and  flow 
measurement  methods  (up  to  20  characters) 

- contians  the  sampling  and  flow  measurement 
times  for  the  methods 

-contains  the  set  up  times  for  the  methods 

- contains  the  code  numbers  of  the  sampling 
methods  (first  method  is  1,  second  is  2,  etc.) 
to  use  at  each  branch 

- contains  the  code  numbers  of  the  flow  measurement 
methods  to  use  at  each  branch  (zero  means  no  flow 
measurement) 


Statements 

I-  7 Subroutine  definition,  COMMON  and  DIMENSION  area 
8- 1 C prints  titles 

11-16  reads  in  and  prints  out  the  sampling  method  names,  times,  and 
set  up  times 

17-18  reads  in  the  method  to  be  used  at  each  branch 

19  DO  loop  to  look  at  all  branches 

20  sets  flag  (used  to  indicate  if  a measurement  is  made  at  a branch) 

21  initializes  maximum  sample  frequency  to  zero 

22  DO  loop  to  look  at  parameters 

23  i f no  measurements  were  allocated  at  branch  I , go  to  15 
(statement  28)  to  check  flagged  points 


I 
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2^-25  selects  the  parameter  number  and  if  its  sample  frequency  is 
greater  than  the  previous  maximum  set  the  maximum  (SMAX)  to 
this  new  frequency 

26-27  set  the  flag  to  indicate  that  a measurement  was  made  at  the 
branch,  and  CONTINUE 

28  if  there  are  not  flagged  points  go  to  18  (statement  ^1) 

29  DO  loop  to  look  at  all  flagged  points 

30-32  obtains  parameter  number  (NSUB) , level  number  (MM) , and  point 
number  (NN) 

33"3^  obtains  the  branch  number  and  if  it  is  not  the  same  as  the 

branch  currently  under  consideration,  go  to  16  (statement  37) 
35”37  set  the  flag  to  indicate  a measurement  was  made  and  determine 
if  the  new  sample  frequency  is  greater  than  the  previous 
maximum.  CONTINUE 

38  i f no  measurements  were  made  at  a branch,  go  to  17  (statement  ^2) 
39-41  obtain  the  sample  method  used  for  the  branch,  calculates  the 

time  required,  and  adds  it  to  the  previous  total 
42-45  print  out  the  appropriate  information 

46  CONTINUE  (look  at  another  branch) 

47"48  pr ints  a title 

49-51  reads  in  and  prints  out  the  flow  measurement  method  names, 

tines,  and  set  up  times 

52  reads  in  the  flow  measurement  code  for  each  branch 

53  DO  loop  to  look  at  all  branches 

54  i f no  flow  measurements  are  made  at  a branch,  go  to  52 

(statement  61 ) 
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S5-57 


i 


1 


r 
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58-60 


6 1 -63 


64-65 

BLOCK  DATA 


obtain  the  proper  location  in  the. SETUP  array  (4  is 
added  because  there  are  4 sample  methods),  calculate 
the  time,  and  add  it  to  the  total 

print  out  flow  measurement  information  and  look  at  the 
next  branch 

print  out  flow  measurement  information  and  look  at  the 
next  branch  (no  flow  measurements) 

RETURN  and  END  control  statements 
- This  segment  is  needed  to  initialize  variables  with 
DATA  statements  that  appear  in  COMMON. 


I 
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SUBROUTINE  PRPAR  - This  subroutine  outputs  parameter  measurement 


f 

I 

[ 


1 

[. 

f 

i: 


L 


I 

I 


t 


information  lising  methods  used,  measure 
locations,  expected  values,  and  equipment  and 
analysts'  times. 


Variable  Definition 
IPPT 
IFPP 

FP  (400,  4) 


POLN  (MPARM,  5) 

I 

AMAR  (MPARM,  MPL) 


NPLA  (25) 

IBN  (51 , MBRPl ) 


- the  total  number  of  flagged  points 

- the  number  of  the  next  flagged  point  to  check 

- contains  flagged  point  information.  FP  (l,  l) 
contains  the  parameter  number  for  flagged  point 
I.  FP  (l,  2)  contains  the  level  number,  FP 

(I,  3)  contains  the  point  number,  FP  (1,4) 
contains  the  method  number 

- contains  the  equipment  names 

- the  current  parameter  number 

- array  containing  selected  methods.  AMAR  (l,  l) 
contains  the  method  number  for  parameter  I. 

AMAR  (l,  2)  contains  the  level  number  at  vi/hich 
the  method  is  used,  and  AMAR  (l,  3"25)  contains 
the  number  of  the  method  at  the  corresponding 
point  in  the  level  is  the  concentration  is 
greater  than  zero. 

- contains  the  number  of  points  at  each  level 

- array  containing  the  branch  numbers  for  each 
level.  IBN  (l,  l)  contains  the  number  of 
points  (branches)  in  level  I.  IBN  (l,  2-50) 
contains  the  corresponding  branch  numbers  For 
each  point  in  the  level. 
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PC  (MPARM,  25,  MBRNC) 


- parameter  concentration  at  a point  and  a 


level  (eg.  PC  (2,  k,  7)  is  the  concentration 
of  parameter  2 at  point  7 of  level  4). 

SAMFRE  (MPARM)  - tne  sample  frequency  of  each  parameter 

PMDATA  (MPARM,  3,  k,  5)  " array  containing  resource  information  for 

up  to  3 methods  per  parameter.  The  information 
is  stored  as  follows:  first  column  - equipment 

codes  (up  to  k items  for  each  method),  second 
column  - equipment  time  per  sample,  third 
column  - analysts'  times,  fifth  column  - cost 
per  sample,  set  up  time,  location  of  analysis, 
and  time  constraint 

SMEQI  I (MEQ.)  - sum  of  equipment  time  for  each  item 

EQTIME  (MEQ)  - total  time  available  for  each  equipment  item 

SUMM  (6)  - contains  totals  of  van  space,  analysts'  times 

and  cost 

Statements 

1-  7 subroutine  definition,  COMMON  and  DIMENSION  area 
8-  I 1 prints  out  titles 

12-13  initializes  LI,  to  indicate  if  there  are  any  allocated  points. 

and  NPTS  to  count  the  total  number,  al located  and  flagged 
1^4-15  finds  the  allocated  level  at  which  the  current  parameter  is 

measured  and  obtains  the  number  of  points  at  that  level 

16  DO  loop  to  look  at  all  points  in  the  level 

17  if  a method  is  not  used  at  a point,  go  to  150  (statement  27) 

18-20  if  it  is  not  the  first  time  through  the  loop,  go  to  100 

(statement  25),  otherwise,  increment  LI  and  NPTS 
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21-24  obtain  the  method  number,  print  out  the  parameter  name, 

method  name,  measure  point,  and  expected  value,  then  look 
at  another  point 

25-27  print  out  the  measure  point,  branch,  and  expected  value  and 
look  at  another  point 

28-29  if  no  points  were  allocated,  just  print  the  parameter  name 
30  initialize  KPP,  the  number  of  flagged  points  for  the  current 

parameter  that  use  a different  method  than  the  allocated 
point  (s) . 

31-32  if  all  the  flagged  points  have  been  checked,  or  if  the  flagged 

point  is  not  for  the  current  parameter,  skip  the  next  segment 
by  going  to  200  (statement  44) 

33-36  obtain  the  level,  point,  and  method  numbers,  and  the  expected 
concentration 

37-39  increment  the  total  number  of  points,  and  print  out  the 

information  about  the  current  point 
40  increment  the  flagged  point  array  location 

41-43  if  the  last  flagged  point  used  the  same  method  as  the 

allocated  point,  don't  increment  KPP.  If  the  method  was 
different,  increment  KPP.  In  either  case,  go  check  the 
next  flagged  point. 

44-46  calculate  the  total  number  of  samples  tal.en,  and  print  it  out 

47-48  print  out  a title 

49  obtain  the  allocated  method  number 

50  DO  loop  to  look  at  all  possible  equipment  items  for  a method 

51-52  if  no  equipment  is  used,  go  to  355  (statement  60) 
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53"55  calculate  the  time  used  for  a method  and  see  if  the  total 
time  allowed  has  been  exceeded 

56-59  print  out  equipment  name,  time,  and  if  a violation  occurred, 
then  look  at  the  next  equipment  item 

60  if  no  flagged  points  used  a different  method,  go  to  kOO 

(statement  60) 

61-62  print  a label 

63  determine  the  method  used  for  the  flagged  point 

6k  DO  loop  to  look  at  all  possible  equipment  items  for  the  method 

65-66  if  no  equipment  is  used,  go  to  400  (statement  75) 

67-69  calculate  the  time  used  for  an  equipment  item  by  a method  and 

determine  if  the  total  time  allowed  has  been  exceeded 

70-72  print  out  equipment  name,  time,  and  if  a violation  occurred, 
then  look  at  the  next  equipment  item 

73-76  print  out  labels 

77-78  obtain  method  numbers  for  allocated  and  flagged  points 

79  DO  loop  to  look  at  each  analyst 

80  calculate  time  used  for  an  analyst 

81-82  if  there  are  no  allocated  points,  go  to  600  (statement  85) 

83  if  an  analyst  is  used  for  an  allocated  method  calculate  the 
t ime 

84  if  the  flagged  and  allocated  methods  are  the  same,  go  to  610 
(statement  86) 

85  calculate  the  analyst  time  for  the  flagged  point 

86-88  determine  if  an  analyst  constraint  was  violated,  and  convert 

minutes  to  hours  and  minutes 

89-91  print  out  analyst  number,  time,  and  if  a violation  occurred, 

and  look  at  the  next  analyst 

92-93  RETURN  and  END  control  statements 
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SUBROUTINE  USECT  - This  subroutine  calculates  the  total  number  of 


items  an  equipment  item  has  been  used. 

Variable  Def i n i t i on 


E (MPARM,  MEQ)  - contains  the  code  numbers  of  each  equipment  item 

that  has  been  used  for  each  parameter 
NP  - number  of  parameters 

NE  - number  of  equipment  items 

U (MEQ.)  - contains  the  number  of  times  each  item  of 

equipment  is  used 

Statemen  ts 

1-  2 subroutine  definition  and  DIMENSION 
3-  4 zeros  cut  array  U 

5-  6 DO  loop  to  check  every  equipment  item  and  every  parameter 

7 if  an  item  of  equipment  is  used  for  a parameter  increment 

the  total  number  of  times  that  item  is  used 
8-10  CONTINUE,  RETURN,  and  END  control  statements 
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SUBROUTINE  TICHAN  - This  subroutine  converts  time  from  minutes  to 


hours  and  minutes. 


Variable 

T 

IHRS 

MIN 


Def i n i t ion 

- time  in  minutes 

- time  in  hours 

- time  in  mi nutes 


Statemen  ts 

! subroutine  definition 

2 converts  from  real  to  integer 

3 converts  minutes  to  hours  (any  remainder  in  the  di 
is  truncated) 

k calculates  the  number  of  minutes  remaining 

5-6  RETURN  and  END  control  statements 


I 
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vision 
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SUBROUTINE  CORRCT  - This  subroutine  accepts  branch  measurement  data 
and  computes  process  model  parameters  from  these  branch  measurements. 
The  output  compares  the  process  model  parameters  to  computed  process 


parameters.  The  subroutine  is  cailed  from  MAIN  if  model  update  or 
verification  is  desired. 

Variable  Definitions 

- contains  the  topology  matrix 

- number  of  sources  in  a particular  study 

- number  of  parameters  in  a particular  study 

- number  of  branches  in  a particular 

- contains  the  parameters  concentrations  for  NS 
sources  and  NP  parameters 

- lists  the  numbers  (from  the  master  list)  of  the 
parameters  in  this  study 

- contains  the  flow  for  each  branch 

- contains  the  mass  parameters  for  each  of  the 
NB  branches 

- contains  the  flow  for  each  source  in  MGD 

- contains  the  mass  parameters  for  each  source 
[C(J,K)*FL0W(J)] 

- parameter  concentrations  for  each  of  NB  branches 
and  NP  parameters 

- parameter  names 

- program  dimension  for  the  number  of  sources 

- program  dimension  for  the  number  of  parameters 

- program  dimension  for  the  number  of  branches 


A(MS0RS,  MBRNC) 

NS 

NP 

NB 

C(MS0RS,  MPARM) 

NPLIST(MPARM) 

EFF(MBRNC) 

X (MBRNC, MPARM) 

FLOW(MSORS) 
P(MS0RS,  MPARM) 

Y(MBRNC,  MPARM) 

P0LN(MPARM,5) 

MSORS 

MPARM 


8^4 


MBRNC 


nfcSijhi,  i5) 


lists  the  branches  for  which  there  are 


6FlIST(25,20) 


iMI2b,25) 

b I a lei.ien  t s 
1-7 

0 “b 

1 u-  I I 
i2-  19 

^ u 

21 

22  129 

-9 
2b 
- 1 


measurement  data 

- For  up  to  25  branches  and  19  parameters  lists 
parameter  numbers  (from  master  list)  for  which 
measurements  are  taken.  Except  BPLIST(IP,1) 
which  lists  the  number  of  parameters  measured. 

- contains  measure  values  for  flow  and  other  para- 
meters for  up  to  25  branches. 


subroutine  definition  and  dimensions 
reads  out  HBiNP  and  NS 
Reads  i n MESUR ( 1 ) 

Reads  in  6PLIST(IB,N) , YM(IB,N) 

Sets  UM  = 0 counts  the  uniquely  deterniined  sources 
Sets  NUN  = 0 counts  the  non-uniquel y lestimatedi 
determined  sources 

Scans  the  branches  for  measurements 
If  no  measurement,  next  branch 

If  all  sources  have  been  determined,  read  ou!  summary 
IFLAG  = 0 treatment  flag 

1C  = 0 counts  sources  contributing  to  branch  it 
measurements  were  taken 

counts  and  identifies  sources  contributing  to  brancn 
if  measurements  were  taken.  If  any  source  flow  has 
been  treated  upon  arriving  at  branch  I, IFLAG  = 10 
It  only  one  source  contributes  to  branch  go  to  I 1 i I ' 
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ko-kl 


'I 


h8 

49 

50-51 

52-56 

57-59 

60 

61-69 

70-71 

72 

73-76 

77-81 

82 

83 

84 

85-91 


Count  the  contributing  sources  which  have  not  been 
previously  determined 
If  only  one  go  to  15(84) 

If  equal  to  zero-next  branch 
If  more  than  one-go  to  next  branch  for  now 
If  the  one  contributing  source  has  already  been 
determined,  go  to  next  branch 

If  not,  set  MESUR(l)  = 2 add  one  to  UN  and  identify 
newly  determined  source  UNIQ(UN) 

Set  IBPLST  equal  to  the  number  of  parameters  measured 
at  branch  1 

Set  flow  of  this  newly  determined  source  equal  to 
the  measured  flow  and  read  out 
identify  parameter  measured  next  BPLIST(|,1P) 
if  treated  go  to  1210(77) 

Set  parameter  of  this  newly  determined  source  equal 
to  the  measured  value  for  branch  1 and  read  out  value- 
if  treated  compute  parameter  of  newly  determined 
source  and  read  out 
look  at  next  parameter  IP 

when  all  parameters  included  in  branch  1 have  been 
considered  consider  next  branch 

arrive  here  if  there  is  more  than  one  source  con- 
tributing to  branch  1 

Identify  the  contributing  source  that  has  not  been 
previously  determined 
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91-93 

94 

95 

96- 

97- 115 


116-117 

118-127 


128 

129 

130 

131- 179 

132- 134 

135 

136-139 

140 

141 

142-143 

144-148 


Add  one  to  UN,  identify  newly  determined  source 
UNIQ(UN) 

Set  MEASUR(I)  = 2 

IBPLST  equal  to  number  of  measurements  on  branch  I 
consider  each  measurement  (IP)  on  branch  I 
if  measurement  is  flow,  subtract  known  contributing 
source  flows  from  measured  flow  in  branch  I and 
read  out  data.  Then  look  at  next  parameter, 
identify  next  parameter 

for  this  parameter,  subtract  known  contributing 
sources  mass  values  from  the  total  parameter  mass 
in  measured  branch.  This  leaves  mass  contributed 
by  new  source.  Convert  to  concentration.  read  out. 
look  at  next  parameter  for  branch  I 
look  at  next  branch 

if  all  sources  have  been  determined  go  to  final 
summat i on 

look  at  each  process  J 

has  the  Jth  process  been  determined  if  so  look  at 
next  process 

set  treatment  flag  I FLAG  = 0 

find  first  branch  in  which  source  J appears 

if  this  branch  was  not  measured  look  at  next  source 

if  branch  was  measured,  is  it  treated,  if  so  set 

I flag  « 10 

set  counters  ICl  » 0 and  1C  » 0 (contributing  sources 
in  branch  I and  branch  (l  - I) 

count  and  identify  sources  contributing  to  branc^  I 
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i 

1 

I 

I 


149-155 

156-158 


1 59-166 

167 

168 

169-171 

172 

173 

174-196 

175 

176 


177-180 

181 

182 

183-184 

185 

186-187 


count  and  identify  sources  contributing  to  branch 
11=1-1 

if  the  difference  in  the  number  of  contributing 
sources  In  branch  I and  II  is  greater  than  1,  look 
at  next  source 

count  and  identify  sources  contributing  to  branch 
I but  not  to  II. 

if  this  is  greater  than  1,  next  source 
if  this  is  not  the  source,  J,  under  current  consid- 
eration, look  at  next  source 

if  this  source  is  the  current  source,  J,  then  set 

MESUR(I)  = 2,  UN  = UN  + 1 and  UNIQCUN)  = J 

are  the  same  number  of  parameters  measured  in  1 

and  II.  if  not  go  to  next  source 

if  yes  set  IBPST  = number  of  parameters  measured 

look  at  each  parameter 

if  parameter  is  not  flow  go  to  750(182) 

if  parameter  if  flow  compute  flow  of  newly  determined 

source  as  the  difference  in  measured  flow  in  branches 

I and  I I . 

read  out  flow  calculation  results 
next  parameter 

if  parameter  is  not  flow  continue  from  here 
identify  parameter  BPLIST(IP) 
if  treated  go  to  790(191) 

if  not,  calculate  the  mass  contribution  of  newlv 
determined  source  as  the  difference  in  mass  from 
measurements  at  branch  I and  I I . Convert  to  con- 


centratIon. 
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188-189 

190 

191-19^ 


195-196 

197 

198 

199 

200-2^9 

201 

202-208 

209 

210-223 

224 

225 

226-246 

111 

228-232 


read  out  results  of  calculations 
next  parameter 

arrive  here  if  treated.  Compute  mass  contribution 

of  new  source  as  the  same  I contributed  by  model 

for  each  parameter  (IP) 

read  out  results,  look  at  next  parameter 

after  all  parameters  look  at  next  source  J.  after 

all  sources  have  been  scanned,  continue. 

if  all  sources  have  been  determined  go  to  summary 

read  out 

set  counter  of  nonunique  (NUN)  or  estimated  sources 
scan  all  branches  again 

if  branch  I was  not  measured  or  if  measurement  was 
used  previously  go  to  next  branch 

otherwise,  count  and  identify  contributing  sources, 

1C  and  ICONS(IC) 

set  counter  NUNPR  = 0 undetermined  branch  sources 
count  and  identify  undetermined  sources  in  current 
branch  and  in  system 

if  no  newly  undetermined  source  go  to  next  branch 

IBPLST  equals  parameters  measured 

look  at  each  parameter  IP 

if  IP  is  not  flow  go  to  561(235) 

if  flow,  compute  flow  of  each  contributing  undetermined 
source  as  PER  times  flow  measured  in  branch  I.  PER  is 
the  decimal  fraction  contributed  by  that  source  in  tf'e 
model,  read  out  results 
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1 


233 

23h 

233-2kk 


2kS 

2A6 

247 

248 

249 

250-260 

261-262 


next  contributing  source 
next  parameter  after  flow 

for  parameter  k,  compute  mass  of  each  contributing 
undetermined  source  as  PER  times  mass  measured  in 
branch  I.  PER  is  the  decimal  fraction  contributed 
by  that  source  in  the  model.  Convert  to  concentra- 
tion and  read  out. 
next  contributing  source 
next  parameter 

total  of  unique  and  nonunique  determined  sources, 
if  all  sources  have  been  determined  go  to  final 
summat i on . 
next  branch 

list  the  uniquely  determined  sources  and  list  the  non- 
uniquely  determined  sources, 
return-end . 
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Flow  Chart  for  Subroutine  CONCK  (feasible  method  determination  only) 
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